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ABSTRACT 

The diversity of Oe icrs pastel system performance 
prediction models currently used by different services 
COntlice With, Ene cOoncepe sot moint operations’ where all 
services must share the common resources to survive. In 
this respect this study presents an analysis and a 
comparison of two operational performance models, the U.S. 
Army’s ACQUIRE and the infrared module of the Navy/Air 
Force Tactical Decision Aid (TDA), WinEOTDA. Differences in 
the modeling of underlying physical principles, input 
parameters, and treatments are analyzed. A comparison of 
the predicted detection ranges is made using a data set 
collected in the Gulf of Oman as the meteorological input. 
Suggestions are sought for the modification of the codes 
that will lead to the same outputs. Finally the possibility 
of adopting one of the codes as a standard TDA is analyzed. 
For the same scenario inputs and with a user-defined sensor 
model WinEOTDA predicted longer ranges for 100% of the 
time. WinEOTDA was observed to be more accurate in 
predicting detection ranges than ACQUIRE because of the 


improved target modeling. 
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toe INTRODUCTION 

Thermal imaging systems are often used for detection, 
BpeCOGMMm= Ton, and dene mhmecmlenw Or targets From ground 
based or aerial platforms by the military. The availability 
of these performance predictions to a decision maker or an 
operator in advance or at the time of operation has vital 
importance for the timely deployment of weapon systems on 
the battlefield. A reliable prediction of performance in 
the target area is also very significant in the mission- 
planning phase of a tactical operation. Tactical Decision 
Pads (TDAs), which can have various forms such as 
nomographs, manuals and computer codes, are tools currently 
used for these purposes to provide predicted detection and 
lock-on ranges to decision makers or operators. The 
performance predictions are currently available in the form 
of computer codes from either the Naval Research Laboratory 
(NRL), the Air Force Research Laboratory (AFRL), or the 
U.S. Army Night Vision and Electronic Systems Directorate 
SIVESD) . 

Current and future modern warfare, which utilizes high 
technology in every means to have the desired impact on the 


enemy, can be analyzed within the concept of joint 


operations. As opposed to the old style battlefield where 
each service had itS own opponent, today every service 
BPequITes Joint resources ana jesme SUPPOre EOnsurVvi Vesna 
requires a cooperative effort, which leads to the concept 
of joint operations. In this respect the existence of two 
different infrared system TDA programs currently used by 
the military conflicts with this idea. This work will seek 
a solution to this problem by comparing the infrared 
modules of the Navy/Air Force TDA, WinEOTDA Version 1.3.3 
dated 1998 and the Army FLIR TDA, ACQUIRE Version 1 dated 
1995, with respect to different means the programs use to 
model target, atmosphere and sensor. 

The objective of this thesis is to determine the 
adifferences in the modeling of underlying physical 
principles, in the input parameters, and in the predicted 
target detection ranges; provide suggestions for 
modification of the codes that will lead to equivalent 
outputs for the same inputs. Finally the possibility of 
using one of them as a standard TDA for all services will 
be examined. We will start Chapter II by presenting some 
fundamentals of infrared radiation theory. This will be 
followed by an analysis and comparison of the ways in which 


this theory is implemented by the two programs. Then the 


—— a 


analysis of the results will be presented in Chapter IV. 


Finally Chapter V will summarize and conclude this work. 


THES PAGE INTENTIONAELEY LEFT BLANK 





II. THEORY AND BACKGROUND 
The purpose of this Chapter is to give a short summary 
of those basic principles of Infrared Radiation which are 


related to the topics addressed in this thesis. 


A. ELECTROMAGNETIC AND IR SPECTRUM 


The electromagnetic spectrum can be described in terms 
of propagating electromagnetic fields that are 
characterized by frequency and amplitude. The optical 
spectrum can be defined as that subset Ot the 
electromagnetic spectrum covering optical wavelengths. 
However there are no exact boundaries for the separation of 
these wavelengths. 

The optical spectrum covers the ultraviolet (UV), 
visible, and infrared (IR) portions of the electromagnetic 
spectrum. Figure 2.1 shows the electromagnetic spectrum and 
identifies various sub-regions of the optical spectrum. It 
can be seen that the visible light spectrum bounds the 
infrared region on the short-wavelength side and the 
microwave bounds it on the long-wavelength side. The 


ultraviolet portion ranges from about Olgas ie) Osa 


micrometer while the visible portion is from approximately 


0.38 to 0.76 micrometer in wavelength. 
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Figure 2.1 - The Electromagnetic Spectrum “From [Ref. 6]”. 


The Sinicaieec Poumon 1S further digmidedmimto moeum different 
sub-regions; the near infrared or short-wavelength infrared 
(SWIR) region (from 0.77 to 3 micrometer), the mid- 
wavelength infrared (MWIR) region (from 3 to 8 micrometer), 
the long-wavelength infrared (LWIR) region (from 8 to 14 
micrometer), and the far and extreme infrared regions (from 
14 to 1000 micrometer) respectively [Ref. 5]. 

Imagers operating in the infrared region of the 
electromagnetic spectrum sense the radiation emanating from 
the targets and the background scene. Unlike night vision 
devices working in the near infrared region, which sense 
the ambient radiation reflected from the targets and the 


backaqmemnd, thermal devices (e.g., Forward Woo kane 


Infrared, FLIR) basically take advantage of the thermal 
energy emitted by the objects in the infrared to detect the 


Signatures. 


B. THERMAL RADIATION LAWS 


Tt is necessary to define some important parameters to 
clarify the basic laws of thermal radiation. The following 
definitions are taken from Seyrafi [Ref. 5]. 

e Absorptivity (@): the ratio of the absorbed radiant 

power to the incident radiant power. 

e Reflectivity (9): the ratio of the reflected 
radiant power to the incident radiant power. 

e TransmisSivity (T): the ratio of the transmitted 
radiant power to the incident radiant power. 

e Emissivity (€): the ratio of the radiant power 
emitted per unit area from a surface to the radiance 
emitted per unit area from a blackbody. 

e Blackbody: defined as an ideal body or surface that 
absorbs all radiant energy incident upon it at any 
wavelength and at any angle of incidence, so that 


none of the radiant energy is reflected or 


transmitted. Blackbodies also have emissivity equal 
to one (€=1). 

e Gray body: a radiation source with an emissivity 
less than unity, and the emissivity is constant over 
all wavelengths [Ref. 1]. 

Table 2.1 gives basic definitions of a few most 


commonly used radiometric quantities. 


Name | Symbol | Units | Description 


| aie sntatned in a red 


contained in a radiation 
fe bee O) 

Radiant ° = ae power traversing 

Flux (Power) a surface. (dQ/dt) 


Radiant flux _ ed. 
Density 
pom “ie 


(Exitance) 
Radiant 


Intensity 
























Radiant flux leaving an 
infinitesimal area of 

surface divided by that 
area. (d0/dA) 


Radiant power per unit 
area incident ona 
(90/94 ) 
Radiant power leaving a 
Point Source per unit 
Solid Angle. (d/dQ ) 
Radiant power leaving or 
arriving at a surface at 
a point in a given 
Girection per unit solid 
angle and per unit area 
projected normal to that 


direction. (07@/dAcos@Q ) 
Table -2.1 = Rag@iemetric Units “Afterei[ Ref si. 
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1. Planck’s Law 


This law gives the spectral distribution of radiant 
emittance of a blackbody radiation source, and can be 
formulated as: 


M(A,T)= Re? peek -] 


where: 
Mijety= the  olkackbody Sspeemeal radian: emittance at 


wavelength > (Watt/cm? um) 


Bay 41 cea O- eWat tlie em 


C1 = 
Co = 1.4388 x 10* um-Kelvin 

a = absolute temperature of the blackbody (K) 
r = wavelength (m) 


2.Wien’s Displacement Law 


Wien’s law 1s simply the derivative of Equation 2.1 
and gives the peak wavelength of the spectral emission for 
a given blackbody temperature by: 

Meow fh ssi) 
where: 


Amax = wavelength where the peak of radiation occurs 


a = temperature (K). 
As the temperature of a source increases, Bias 
equation indicates a shift in the wavelength of the maximum 


radiation toward a shorter wavelength. This can be observed 


Gaspmbedlivy in Figure 222. 
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Figure 2.2 - Spectral Radiant Exitance of Blackbodies at 
Various Temperatures “From [Ref. 5]”. 


3. Stefan-Boltzmann Law 


This law is simply the integral of Equation 2.1 and 
provides the total radiant emittance by integrating 
Planck’s law over the entire spectrum. The following 
equation applies only to blackbody and graybody sources 


[Ref. 1], and is formulated as: 


M =eoT* (2235 


1 








where: 
€ = emissivity 


M 


total radiant emittance of a blackbody 


oO the Stefan-Boltzmann constant (5.7 x 10° Watt/m’*) 
4. Total Power Law 
The radiation incident upon a body may be transmitted, 
absorbed, or reflected and by conservation of energy the 


sum of the ratios of each of these components to the 


incident power must be one. 


A+ p+t=l1 | (2.4) 
where: 
@ =Sabsonmp ta wus, 
p = reflectivity 
Te= tCancmussivaty 


5. Kirchoff’s Law 
This law states that the bodies emit as well as they 


absorb at any wavelength; this can be expressed as [Ref. 


Z|: 


JaIL 


Ce ATMOSPHERIC PROPAGATION 


The optical signal radiated from an object must pass 
through the Earth’s atmosphere before it reaches a 
receiver. No matter how strong the target Signature is, the 
intervening atmosphere always attenuates the thermal 
Signal. This attenuation 1s due to the individual or 


collective effects of the following phenomena: 
e Molecular absorption, 
e Molecular scattering, 
e Aerosol absorption and, 


e Aerosol scattering. 

Molecular absorption is due to the ability of certain 
molecules to go from one vibration-rotation state to 
another. ite 1s generally characterized by discrete 
absorption lines arising from the quantal nature of the 
absorption, modified by broadening processes, including 
pressure and Doppler broadening. Water vapor is the most 
important of these molecules. It limits the useful range of 
infrared wavelength to the 3-5 and 8-14 micrometer bands. 
Other molecular absorbers include carbon dioxide, carbon 


monoxide, ozone, methane and nitrous oxides. 


sa 





Scattering, on the other hand, is the redistribution 
of the incident energy into all propagation directions. The 
scattering processes are generally related to the ratio of 
wavelength to particle size, resulting in broad absorption 
spectra, generally maximum when wavelength 1s pi times 
particle size. When wavelength is much shorter than 
particle radius, Rayleigh scattering is dominant; when 
wavelength is much longer than particle radius scattering 
becomes small. Scattering differs from absorption in that 
the scattered radiant energy remains in the same form as 
the incident radiation [Ref. 1]. Water droplets suspended 
in the air are the most important source of scattering. 
Other sources of scattering include dust, smoke, smog, 
isc, Or SNOW [Reft. 10]. 

Eerie On = Ei = sine "OF AalSsorpe tommeanc Scatter inea: == 4S 
the process of attenuation of the radiant flux in passing 
through the atmosphere. It can be expressed in terms of an 
exponential coefficient used in the following formula, 
called Beer’s law, where p is the extinction coefficient 
and R is the path length: 

Te 
T is the transmittance of a path length R through the 


atmosphere. The spectral transmittance for 1 km path length 


ALS 


at sea level under “typical” conditions 1s shown in Figure 
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Figure 2.3 - Typical Atmospheric Transmission for a 1 Km Path 


Length “Fronmether. 1 ].” . 
It can be observed that atmospheric extinction 1s a strong 
function of wavelength, which severely affects the 
transmittance through the atmosphere, and 3-5 wm (MWIR) and 


8-12 um (LWIR) wavelength regions are the only ranges valid 


for atmospheric propagation. 

The large number of parameters involved in optical 
transmission Enrough the atmosphere makes numerical 
calculations of atmospheric transmission inevitable. The 


aim of the numerical calculation is to predict with a high 


degree Oye accuracy the transmittance iE aa @ Uicia the 
atmosphere, given a path, atmospheric Eondilt2ons- 
wavelength, and a set or measured Or predicted 
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meteorological parameters (Ref. Srllee ene an accurate 
atmospheric transmission cCalculaeon all molecular, 
aerosol, and precipitation effects must be considered, and 
a detailed model must be used to get precise results when 
necessary. Three such models are LOWTRAN [Ref. 17], FASCODE 
[Ref. 3], and MODTRAN [Ref. 18] that are used to obtain the 
atmospheric transmittance T(A). 

The SEARAD Radiance Model, which is a surface radiance 
model integrated with the MODTRAN2 transmission model, will 
be used to predict the atmospheric transmittance required 
by the sensor performance model in this thesis, and will be 


described in that context. 


D. TARGET SIGNATURE 


Infrared sensors respond to the difference in radiance 
between target and background. From Equation 2.3 it is 
evident that the target-background radiance difference can 
be related to an equivalent temperature difference (AT), 
which appears as a thermal quantity. The equivalent 
temperature difference is defined as the temperature 
difference of two blackbody sources required to produce the 


actual radiance difference between target and background. 


ils 


Figure 2.4 shows the geometry, which can be used to obtain 


AT for an extended target. 


Extended source 
target 


Area of target Infrared sensor 


Radiance seen by the detector ” Colle | eee 








or eee eee ee et ee eee 


! ze 
| Range, R 
Figure 2.4 - Differential temperature geometry “From 
[Ref.1]”. 


However, a more important quantity than temperature 


yc 
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This is the equivalent blackbody temperature difference 
seen through an atmospheric path that produces the same 
sensor output voltage difference as the real target and 
background. The following figure pictures the difference 


between the concept of temperature difference at zero range 


wave 
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aperture of the sensor (apparent delta T, OMT a jee 
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Figure 2.5 - Apparent delta T “From [{Ref. 1]”. 


There are various computational techniques available 
to determine the apparent target-to-background temperature 
Gifference at the entrance aperture of a broadband infrared 
sensor. In this thesis the following two techniques will be 
used for calculating the apparent delta T where necessary. 

The first technique estimates an apparent temperature 


difference as the product of the target-to-background 


temperature CAL. and the atmospheric broadband Beer’s law 


Eransmittance (@,), Using an extinction coefficient 


averaged over the system bandwidth. That is: 
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The broadband transmittance (7,,) is defined for a one- 


kilometer path length and R is the target-to-sensor range 
in km. 

While Beer’s law is valid for monochromatic (single 
wavelength) sources, infrared imaging sensors typically 
operate with a broad bandwidth of several micrometers. 
Propagation of broadband radiation presents significant 
computational difficulty, since Beer’s Law is not generally 
valid for broadband transmission of light [Ref. 1]. In the 
broadband Beer’s Law approximation, a band averaged 
extinction coefficient is computed from the transmittance 
at a reference path length. The transmittance is found by 
averaging the spectral transmittance over the wave band for 
that path length. In this computation the reference 
extinction coefficient is then taken to be constant over 
that bandwidth for all ranges. However, in actuality the 
spectral extinction coefficient varies within the bandpass, 
and the band averaged extinction coefficient will be a 
function Of the Yrange. Thus an broadband transmission, 
absorption is not characterized by a constant exponential 
coefficient as in Beer’s Law, and the exponential! mange 


dependence does not hold. 
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The second technique™is differeénetrom *the Viivrst ih 
that a broadband Beer’s law assumption is not used to 
determine the atmospheric transmittance. Instead an 
atmospheric transmission program is used to find the 
broadband transmittance directly as a function of range. 


Then as in the following equation the product of target-to- 


Pe itGeemm@eamecemperazure (Ale and eehe output transmittance 


values of the atmospheric transmission program, a(R), is 


taken to determine the apparent temperature. 


Mae cayiNl_ aalle) (2.8) 


It must be noted that differential target temperature (AT...) 


tet 
used in Equation 2.7 and 2.8 is referenced to two blackbody 
sources required to provide the same differential flux as 
that of actual target and background. Thus the temperature 
is not in fact attenuated through the atmosphere; energy or 
radiance is attenuated [Ref. 1]. 

This technique is the one presently used in the system 
performance program ACQUIRE [Ref. 16] in this thesis. This 
technique eliminates the errors associated with a broadband 


Beer’s law assumption [Ref. 1]. 
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E. DETECTION CRITERION 


Target detection refers to different levels ene 
distimngqurcshing an™@object from background: The lowest levem 
1s simply a detection of the object. The highest level is 
the identification of a specific object. These levels can 
be gathered sigh aqs, two groups: pure detection and 
discrimination detection. In pure detection locating an 
object in the scene is sufficient to declare detection. On 
the other hand, in discrimination detection where the scene 
contains many non-targets, objects cannot be detected as 
targets until sufficient shape information can be obtained 
to distinguish the target from non-targets or clutter. 

The traditional FLIR analysis describes the 
interaction of the FLIR-aided eye with two types of simple 
targets: an isolated rectangle, characterized by uniform 
temperature difference from the background and a periodic 
bar pattern, also characterized by a temperature difference 
from the background [Ref. 10]. The minimum temperature 
difference required for detection of the rectangle is known 
as the minimum detectable temperature difference (MDTD). 
The temperature required to resolve the four bars is known 


as the minimum resolvable temperature difference (MRTD). It 


8 





is common to represent operational targets with equivalent 
bar targets for evaluation. Johnson [Ref. 19] conducted a 
number of experiments with a number of trained observers at 
the U.S. Army’s Night Vision and Electronic Sensor’s 
Directorate to develop resolution requirements Mee 
detection, recognition and identification of objects by 
using these bar patterns. He determined the average number 
of line pairs required for different discrimination levels 
as listed in Table 2.2. Today these are known as the 
myehnson “Criteria”. The cycle ~Griteria in Table 2.2 
correspond to the necessary number of resolution elements 
on the critical dimension of the object with a  two- 
dimensional cycle requirement and to a probability of 50% 


for a Given discrimination task. In this tablew=n, is 


the number of cycles required to be resolved in order to 


achieve a 50% probability of discrimination. 


Detection An object within the sensor FOV is a target 
Ns, =0.75 of potential military interest 


Classification The target belongs to a general class of 
Ns =1.5 vehicles: tracked or wheeled 






Recognition The target is a specific object within a 
Nsy =3 class of similar objects: tank or APC 


Tdentification The target is a specific vehicle:T72 
Nsy =6 | 


Table 2.2 — Johnson Cycle Criteria “From [Ref. 13]”. 
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In the detection process the Johnson criterion is used 
after finding the target, where the target size and shape 
provide information Oe detection, mecognition, and 
identification. It provides the connection between the MRTD 
and field performance of the sensor. 

In two-dimensional discrimination, also used in this 
thesis, target area is more important than the minimum 
dimension used in one-dimensional detection, as first used 
by Johnson. The “critical dimension” as used in two- 
dimensional resolution is defined as the square root of the 
Carget area. 

The two-dimensional FLIR92 model uses the critical 
dimension approach in the same manner [Ref. 9]. In this 
thesis Shumaker’s [Ref. 10] approach, which takes into 
account the aspect angles will be used to calculate the 


critical dimension of the target. 


F. FLIR SYSTEM PERFORMANCE PARAMETERS 


The Eolmbow ing section describes the phys 1Gam 
parameters which determine the MRT and MDT. 


1. Physical Parameters 


The physical parameters of a FLIR system are defined 


as. follews: 


ae 


a) Field-of-View (FOV) 
The FOV of an IR system is the angular space in 


which the system accepts radiation. The system FOV and the 
distance, or range, from sensor to the object determine the 
area that a system will image [Ref. 1]. 

b) Instantaneous Field-of-View (IFOV) 

The instantaneous FOV is the angular cone through 
which a detector senses radiation; i1t depends upon the 
optical design. It includes both the optical blur diameter 
and the DAS. When the blur diameter is small compared to 
the DAS, the IFOV and DAS are approximately equal [Ref. 9]. 

c) Detector Angular Subtense (DAS) 

The detector angular subtense is used to describe 
the resolution limitations of the detector size. DASs in 
the in-scan (Ax) and cross-scan (Ay) directions are given 
by the detector width or height divided by the focal 
length. It describes the best resolution that can be 
achieved by an EO system due to the detector size 
intimitations [Ref. I]. 

d) Modulation Transfer Function (MTF) 

The MTF of a system iS a primary measure of the 


overall system resolution. The system MTF gives’7 the 
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transfer of input spatial frequencies, and it can be in 
both horizontal and vertical directions. 

The modulation transfer function is the magnitude of 
the optical transfer function, which actually alters the 
image as it passes through the optics and circuitry of the 
system. It can be formulated as the output modulation 
produced by the system divided by the input modulation at 


that spatial frequency: 


MTF = OUTPUT _MODULATION 
INPUT __MODULATION 


(2 diGy 

2. Noise Equivalent Temperature Difference (NETD) 

NETD is the temperature difference between a large 
target and its background, which produces a SNR of one in 
the video signal. In performance predictions, it 1S used as 
an intermediate sensitivity parameter for simplification of 
formulations of performance parameters such as MRT, and 
MDT. NETD can also be described as a system’s ability to 
detect small signals in noise. It does not account for the 
Spatial and temporal integration effects of the eye. 


Shumaker [Ref. 10] gives NETD for a scanning system 


aS: 
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10(FOV, FOV,F.N,.N,,}? 
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where: 

D is the aperture diameter (m) 

De is the band average detectivity of the detector 
with no cold shield (cm Hz” W’*) 

Np is the number of detectors 

Nsc 1s the scan efficiency 

Ax 1s the in-scan detector angular subtense (mRad) 

Ay is the cross-scan detector angular subtense 
(mRad) 


ON/dTis the derivative of Planck’s Law (the “Thermal 
Meacwanie’ — Watters 6k sr) 

FOV; 1s the in-scan field of view (mRad) 

FOV, is the cross-scan field of view (mRad) 

Nas is the overscan ratio 

Nee is the serial scan ratio 

le is the frame rate 

However the concept of three-dimensional noise, which 
has been successfully integrated into the U.S. Army’s Night 
Vision and Electronics Sensor Directorate’s FLIR92 sensor 


model, will be used for defining the infrared system noise 


nS 


in this thesis. This method eliminates the limitation of 
NETD on defining only the temporal detector noise, and 
characterizes the noise both spatially and temporally from 
various sources. 

3. Minimum Resolvable Temperature Difference (MRTD) 

Minimum resolvable temperature difference (MRTD) is 
the most used and useful FLIR specification parameter. It 
is defined as the temperature difference between the 
background and a set of four standard bars (7:1 aspect 
ratio) required to make the bars just resolvable, as a 
function of the spatial frequency of the bars. [Ref. 10] 

There are several important features of MRTD. First, 
atte 1s an end-to-end system measure Ineludingeboimm 
resolution and sensitivity, and it is subjective since it 
involves the judgment of the human observer. Second, the 
temperature difference that is required to resolve the four 
bars increases as the bars become smaller, as can be seen 
from Figure 2.6. Finally the MRT curve is asymptotic at a 


spatial frequency near 1/DAS, where the MTF becomes zero. 
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cyc/mrad 


Figure 2.6 - MRTD Patterns of Differing Spatial Frequency (&) 
[Ref. 1] 

There are several formulations used and various 

authors have proposed some alternative expressions for the 


MRTD. Shumaker [Ref. 8] gives the following formula: 


20SNRT(FOV, FOV,v?) p? 


0 ss ns SSS 
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where: 


SNRT is the perceived signal-to-noise threshold 


V 1s the spatial frequency (cycles/mRad) 


D 1s the aperture diameter (m) 


PLE | 


Dp* * Sas 


the band-average detectivity of the detector 


with no cold shreVvawiem Hz” w+) 


Np is 
Nsc 18 
Nes 18 
ae is 
DSc 1s 
Ay 1s 
(mRad) 
MTFs is 
L 1s 
tee is 
ON/OT is 


(or. is 
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the 
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number of detectors 

scan efficiency 

cold shield efficiency 

transmittance of the opti¢s 

in-scan detector angular subtense (mRad) 


cross-scan detector angular subtense 


system modulation transfer function 
length-to-width ratio for the bar chart 
eye integration time (0.2s) 


thermal derivative of Planck’s Law (watt 


in-scan field of view (mRad) 


cross-scan field of view (mRad) 


mouse Lilter facror 


4. Minimum Detectable Temperature Difference (MDTD) 


The description is almost the same as the MRTD of 


FLIR system. 


representation 


rather than a four-bar target. 


a 


The difference between the two 1s the 
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of the target, which for MDTD a _ square 


MDTD of a FLIR gives the 


temperature difference between an isolated square and a 
uniform background that renders the square just detectable, 
as a function of the dimension of square in spatial 
frequency. As in the case of MRT it has an element of 
subjectivity since the judgment of an observer is involved 
in the process. In the observation process the observer 
approximately knows the target location. MDT is given as: 
20SNRT (FOV, FOV, AxAy}7(Q, +72) 


MDT (Vv) >—_$§$££____—__ 23 ) 
mt,DAxAyD™ dN/dT7,,(N p71.) Qr(r2 +72 +2] 


SC t, 


where: 


SNRT 1s the perceived signal-to-noise threshold 

Vv 1s the spatial frequency (cycles/mRad) 

Qe is the solid angular subtense of the target 
(mRad) * 

D is the aperture diameter (m) 

D** ais the band average detectivity of the detector 
with no cold shield (cm Hz” WwW?) 


Np 1s the number or Gerectore 
Nsc 1s the scan efficiency 


Nes 1S the cold shield efficiency 


Zo 


a fs tne transmittance ef ene Cpteres 

Ax °1s the in-scan detector angular subtense (mRad) 

Ay 1s the cross-scan detector angular subtense 
(mRad) 

eg is the eye integration time (0.2s) 

ON/dOT is the thermal derivative of the Planck’s Law 
‘sieikic cio 2 = ete =) 

FOV, is the in-scan field of view (mRad) 

FOV, is the cross-scan field of view (mRadqd) 

Xs 1s the resolution of the system that includes the 
front-end resolution and back-end resolution (mRad) 

ie is the resolution of the back-end that includes 
the detector electronics ee eee Beso tome 
resolution of the multiplexer, resolution of the display, 
resolution of the eye, and the resolution due to image 
motion (mRad) 

As seen from the above equation MDT has no first-order 


linear dependence on MTF, which means that it does not show 


the asymptotic behavior that MRT does. 
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IIIT. MODELS 
i. GENERAL INFORMATION ABOUT TDAS 

Tactical Decision Aids (TDAs) are tools that assist a 
decision maker or an operator in planning or performing a 
task. They can be in such various forms as nomographs, 
manuals or computer codes, which is the form used in this 
thesis. They are designed to aid a decision maker by 
assimilation and convenient presentation oof data and 
analysis of a tactical problem beyond what is feasible by 
humans in timely fashion [Ref. 8]. 

In parallel with the rapid development in technology 
of new weapon systems, it 1s becoming more complex to plan 
or decide on the timely deployment of these systems on the 
battlefield. In order to have the desired impact on the 
targets, TDA codes used by the personnel must be quick and 
user-friendly to accelerate the planning or operational 
process. These models can also be used in the design or 
testing phase of a new system. 

Fach code contains the following three fundamental 
parts; a) Target Model which determines the inherent signal 
emanating from the target and background and converts the 


radiance difference between them into a temperature 


ot 


difference (A7) at zero range. b) Atmospheric Model, which 
1s the module that calculates the apparent delta T by 
estimating the degradation of signal due to the atmosphere 
at the entrance aperture of the sensor. c) Sensor Model 
which describes the sensor woeroonenee im terms of MRITD em 
MDTD as a function of spatial frequency. This model 
determines the detection or the lock-on range of an 
electro-optical system when applied to the apparent target 
Signature. In this thesis the Johnson criterion will be 
applied mene specifying a detection decision where 
necessary. 

The following two sections will cover the models used 
to design and calculate the performance parameters (1.e., 
MRTD and MDTD) of a sensor and the calculation of 
atmospheric transmission, as required by the WinEOTDA and 
ACQUIRE models. In the remaining sections TDAs under study 
will be described according to the fundamental parts listed 


above. 


B. FLIR92 MODEL 


FLIR92 is a system evaluation tool that uses basic 
sensor parameters to predict overall system performance for 


thermal imaging systems. It 1s a desktop computer model 


a2 


working in the DOS environment. The model calculates 


modulation transfer function (MTF), noise equivalent 
temperature difference (NEED) ; minimum resolvable 
temperature difference (MRTD) , and minimum detectable 
temperature difference (MDTD) by using basic system 


parameters. The principal function of the model is to 
predict whether or not a system achieves the required MTF, 
system noise, MRTD, and MDTD determined necessary to meet a 
target acquisition and discrimination task. 

FLIR92 models parallel scan, serial scan, and staring 
thermal imagers operating in the mid and long-wave infrared 
regions. It can be used for thermal imagers only and cannot 
predict the performance of other kinds of electro-optical 
sensors. The model does not predict target acquisition and 
discrimination range performance. [Ref. 11] 

In FLIR92, there are two different outputs: an MRTD 
Gemmonii~e used» fom whaeh a sdascriminatiom decision is emade, 
and an MDTD commonly used for which an acquisition decision 
1S made. 

FLIR92 calculates the system’s overall MTF by using 
linear filter theory. The MTFs for the components are 
multiplied together. Instead of including an MTF for each 


component, MTFs of common system components are gathered in 
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three main Gisoups:: a) Prefilter MTFs, b) Temporal 
Postfilter MTFs, and cc) Spatial Post#iilter M@PFs. The 
components may vary according to the design and the users 
can add new MTFs into these groups. 

FLIR92 ignores signal and noise aliasing in the MRTD 
and MDTD predictions, for thermal imaging systems are 
assumed to be well designed, so that image artifacts due to 
under-sampling do not significantly degrade the system. 
Thus, the model is implemented with enough flexibility to 
accommodate most system designs through user determined 
pre- and post-sampling MTFs. Also, MRTD is not predicted at 
Spatial frequencies exceeding the Nyquist frequency [Ref. 
ie 

As opposed to the first generation thermal imaging 
systems where NETD was used to predict the system 
performance, in second generation systems noise was defined 
in a three dimensional coordinate system (temporal, 
horizontal spatial, and vertical spatial) by the FLIR92 
model. The model calculates the full temporal noise, and 
the spatial noise is incorporated into the MRTD prediction 
via the three dimensional noise model and summary noise 


Seaecors. “NOLSeG  “CcCaleumations “are Made —=relative to a 
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measuring port that is assumed to be located at a video 
port prior to the system display.” [Ref. 11]. 

FLIR92 calculates the horizontal and vertical MRTDs 
depending on the direction of the standard four bar 
pattern. As mentioned in Chapter I the MRTD depends 
directly on the system transfer function, which is 
represented by the system overall MTF, and the system 
sensitivity that is described by NETD. To predict MRTD and 
MDTD, the spatial integration of the eye/brain system must 
be modeled. FLIR92 uses a synchronous integrator model for 
MRTD predictions as opposed to a matched filter model. 
“With this method, the eye/brain system 1S assumed to 
Spatially integrate over the image of a bar, ignoring 
blurring of the target caused by finite apertures in the 
system. ” (Ref. 11). In the case of periodic targets, 
synchronous integrator and matched filter methods give the 
same results. However, since the algorithms required to 
implement are simpler than matched filter algorithms, 
FLIR92 uses the synchronous integrator method. On the other 
hand, MDTD prediction is based on the matched filter 
concept, in which the eye/brain filter is matched to the 


Signal in order to maximize the signal to noise ratio. 
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Ce SEARAD RADIANCE MODEL 


SeaRad is a FORTRAN computer code used in predicting 
the radiance of the ocean surface. It includes a more 
accurate description of the sea surface including effects 
of solar heating and reflection and wind modification of 
the sea surface. For transmission it uses a modified 
version of the U.S. Air Force program MODTRAN2, which uses 
a card input system to compute atmospheric transmittance 
and path radiance. SeaRad is DOS-compatible and runs on a 
personal computer. In this thesis a Matlab shell for input 
and output of this code [Ref. 4] was used to compute the 
atmospheric transmittance values required by the ACQUIRE 
model . 

The SeaRad surface state model is based on the Cox- 
Munk statistical model for wind-driven capillary wave 
facets. It operates exactly like the original MODTRAN2 code 
with an additional new logical parameter that 1s required 
in the aunput file. “Sun glint is “included in the “Sea 
radiance prediction provided that the user has chosen to 
execute SeaRad in radiance mode with solar scattered 
radiance included.” [Ref. 7]. The program is valid for the 


spectral range from the visible to far infrared regions. 
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D. WINEOTDA MODEL 


Windows Electro-Optical Taewi Gal Decision Aid 
(WinEOTDA) iS a computer model developed by the Naval 
Research Laboratory (NRL), in Monterey, CA [Ref. 22]. It 
was derived from the Electro-Optical Tactical Decision Aid 
Mark III (EOTDA III), which was originally an EOTDA model 
running in DOS. The USAF Philips Laboratory first created 
Mark III, and then NRL incorporated Navy sensors into this 
model. NRL developed the Windows version of this program to 
make it user friendly and simplify the prediction process. 


WinEOTDA predicts the performance of electro-optical 


weapon systems and night vision goggles (NVG), working in 
the infrared (8-12 micrometer), visible (0.4-0.9 
Micrometer), and laser (1.06 micrometer) wavelengths region 


of the optical spectrum. The prediction is based on 
environmental and tactical information, which includes 
meteorological data, time over target, target location and 
characteristics, sensor specifications and height, and 
background characteristics. 

WinEOTDA uses a Graphical User Interface (GUI) design 
to display the maximum information on the screen and 


present the inputs and outputs in a single window. This 


OF, 


allows the user to reach the details just by clicking on 
the links on the main screen. Using the main menu and 
toolbar on the same screen can also make the selection. All 
meteorological and operational data can be input via the 
drop down menus and links. 

WinEOTDA consists of three basic components: target 
model, atmospheric transmittance model, and sensor 
performance model. In the following three sections each of 
these components will be explained and a summary 
description of output files will be given. 


1. Target Model 


Target model calculates the strength of the electro- 
optical signal at see range using target and background 
characteristics entered by the user. The radiance 
difference between the target and background 1s converted 
to an equivalent blackbody temperature difference via the 
thermal model Target Contrast Model #2 (TCM2). TCM2 is a 
very powerful and accurate target Signature model developed 
by Georgia Tech Research Institute. It is based on heat 
transfer and treats the target as a distinctive three- 
dimensional network of nodes that exchange heat with one 
another as well as with their environment [Ref. 12]. The 


model provides a very detailed target signature. 
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The TCM2 model ceontinweusly calculates a new 
temperature for all nodes during the period between the 
beginning and the end of the operation at various time 
intervals determined by the user. For each interval TCM2 
computes a mean temperature and identifies the hottest and 
coldest of the visible nodes [Ref. 12]. Then the one with 
the greater contrast to the background is identified and 
the sensor model uses the temperature and projected area of 
this facet in MDTD based range detection. The target mean 
temperature and total projected area are used to compute 
MRTD based detection range. 

The WinEOTDA version 1.3.3 dated 1998, the version 
used in this thesis, includes 20 different targets in its 
target menu, containing land vehicles and buildings, 
aircraft, and naval ships. Target heading, operating state, 
and speed of these targets provide the necessary input for 
the TCM2 model to calculate internal heat sources as well 
as surface heating and cooling. Target heading affects the 
perceptible solar heating on target, as the target 
operating state gives Infermation about the heat 
interaction with the environment, and the surface heating 


of the target. The movement of the target, represented by 


SS, 


1ts speed and the wind speed, provides a cooling effect on 
the target. 

The backgrounds in the model are grouped under two 
different categories, as general backgrounds and specific 
backgrounds. The general background offers five 
subcategories: continental, urban, desert, ocean, and snow. 
It describes the dominant terrain feature of the target 
area, which gives the information used to calculate the 
solar reflection by the model. The background, which is the 
immediate area surrounding the target, consists of eight 
different structures: vegetation, soil, snow, water, 
concrete, asphalt, swamp and rocky field, which are further 
described by the composition, coverage or depth of the 
surface type. Three different backgrounds used in this 
thesis are water, soil, and vegetation, so as to represent 
a beach scenario for joint operations. Despite the use of 
multiple backgrounds in the model, the program uses the one 
entered first as the primary background to calculate the 
solar heating and reflection of the ground. However these 
backgrounds are not considered to be independent and the 
program directs the user to enter the most representative 
one first. In the case of water background the depth 


affects the heat capacity of the water body and clarity 
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affects the heat flow from the surface [Ref. 12]. Soil 
types also define the heat capacity and reflectivity, as 
the moisture affects cooling rate, which is also an issue 
for the vegetation model. 


2. Atmospheric Model 


The atmospheric model calculates the degradation of 
the signal in transit from target to sensor. A limited 
version of the LOWTRAN atmospheric propagation model is 
used in calculations to predict the transmittance through 
the atmosphere. The path radiance is not included in this 
modeling. The model for a range of four kilometers is used 
to evaluate the transmittance and thus the band averaged 
extinction coefficient. Then the Beer’s ene approximation 
1s used to calculate the other transmittance values for 
different ranges. 

WinEOTDA uses the two-layer model, which calculates 
two extinction coefficients for below and above the 
boundary layer height. A weighted average of transmission 
is used for sensors above the boundary layer. 

The aerosol modeling consists of 19 aerosols from 
MOWTRAN 5,6,7 [Ref. LT ie and the Navy Aerosol Size 
Distribution Model (NAM). WinEOTDA aerosol models include: 


ialjtica |. urban, maritime, tropospheric, desert, navy 
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maritime, advective fog, radiative fog, and camouflage 
smokes. The navy maritime model, which is used in this 
thesis, describes aerosols found in the boundary layer of 
oceanic environments. WinEOTDA includes nine different 
aerosols in the Navy maritime model, defined by air mass 
history, and the 24 hour average and local wind speeds, a 
distinguishing factor from the standard maritime model in 
LOWTRAN. 

Meteorological data are input by using the Met input 
screen for the transmission calculations in WinEOTDA. The 
following parameters of the target scene are required by 
the model on an hourly basis: surface temperature, surface 
dew (Oe Gaye temperature, aerosol, battlefield induced 
COMtaminants “—(BlCs) P= visibi wey, “paeelpltallon “index, “sain 
rate, wind speed and direction, boundary layer height, low, 
middle and high cloud data. 

The surface temperature and dew point are used to 
compute the relative humidity. Then relative humidity, 
aerosol and visibility parameters are used to calculate an 
exXternmction cociiverents In Navy Maritime, the model 


calculates visibility. 
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3. Sensor Model 


The sensor performance model evaluates the range at 
which the signal received by the sensor equals’ the 
threshold value for detection. The target apparent size 
(angular subtense) as viewed from the sensor determines 
this threshold value as where the angular subtense is equal 
to the critical dimension of the target divided by the 
range to sensor. 

WinEOTDA supplies the user with a number of sensor 
data files identified by a unique three-digit index. The 
operator selects the sensor according to this number from 
the sensor list. The program offers two kinds of sensor 
IDs: standard IDs reserved for sensors supplied with the 
program and additional IDs for user-defined sensors. The 
physical and performance parameters of the supplied sensors 
are encrypted into separate data files and kept in the 
program in pure ASCII code. The identifications of these 
sensors are not available to the user. In this thesis a 
user-defined model using the standard SADA II scanning 
focal plane array was designed as a second generation FLIR 
sensor by the FLIR92 model using the basic’ sensor 


parameters found in the literature. 
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4. Output Files 

WinEOTDA outputs are displayed in three different 
formats that can be selected via the main menu or toolbar. 
Alphanumeric, graphic, and tabular outputs are created 
automatically after each run. An alphanumeric output is 
designed with the following parameters: MRT Range, MDT 
Range, Lock-on Range, MRT Delta T, MDT Delta T, Lock-on 
Delta T, Background Temperature, MRT Target Temperature, 
MDT Target Temperature and Lock-on Target Temperature. 
Graphic Ouse DWE includes the SUEDUtE range, target 
temperatures, and delta T values while the tabular output 
displays only the output ranges according to different 
viewing directions. In addition to these outputs the model 
displays the maximum ranges for each target on the main 
screen. The units of the ranges can be changed via the main 
screen, which has the options oof kilo-feet came) 
kilometers (km), and nautical miles (nm). WinEOTDA gives 
maximum range predictions for detection only, for both 


Narrow and Wide FOVs. 


E. ACQUIRE MODEL 


ACQUIRE is a range performance program that was 


developed by the US Army CECOM, Night Vision and Electronic 
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Sensors Directorate (NVESD). The version used in this 
thesis is dated May 1995. It runs on IBM compatible 
personal computers in a DOS environment and on Unix 
workstations. 

ACQUIRE predicts target detection and discrimination 
range performance for systems working in the visible, near 
infrared, and infrared spectral bands. There are two 
different range prediction tasks in ACQUIRE: target 
discrimination and target spot detection. Two-dimensional 
Johnson cycle criteria along with MRTD predict the target 
discrimination ranges, while target spot detection (star 
detection) ranges are predicted by utilizing SNR theory and 
moanie MOD. 

As explained in ACQUIRE’s User’s Guide “ACQUIRE is 
intended ven experienced systems analysts who are 
knowledgeable of 1) the principles of imaging electro- 
optical systems and their application to target acquisition 
tasks, 2) the parameterization of target acquisition 
scenarios for the purpose of evaluating targeting systems, 
and 3) the basic methodologies applied in the model.” [Ref. 


eS 
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1. Target Model 

ACQUIRE uses the target information given by target 
Signature and dimensions to calculate the probability of 
target discrimination. The target signature is represented 
by target contrast in the visible or near infrared. regions 
while it is defined to be the temperature difference 
between target and background at zero range in the infrared 
spectral bands. The probability of target discrimination is 
a function of the number of equivalent cycles resolved on 
the target by the sensor. The number of cycles resolved is 
determined by minimum resolvable temperature difference 
(MRTD) for a target at a given range and apparent Signature 


[Ref. 13] and is given by: 


j 
n=—f, (ely 
e 
where 
Ge is the characteristic size of the target (m), 
i 1s the range to the target (km), and 


Fi is the frequency (cycles/mRad) resolved by the 


sensor for the target at range r. 
Then the probability of discrimination is calculated 
by utilizing the following target transfer probability 


function (TTPF) used by ACQUIRE as a curve fit: 
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ft (1/155) Toe 
1+ (n/n, )° 


where 


E is equal to 2.7+0.7(n/ no ) , 


Io 1s the number of cycles required to be resolved 


in order to achieve a 50% probability of discrimination. 

ACQUIRE offers 23 land targets in its target look-up 
table file. For targets that are not represented in the 
internal lookup table, a data file is built with the target 
dimensions and signature. When the target signature is not 
entered in the model, a thermal default of 1.25 degrees C 
moe Used Tomeit. 


2. Atmospheric Model 


ACQUIRE offers two different methods of modeling 
atmospheric transmittance. The first uses the Beer’s Law 
approximation calculated from the atmospheric transmittance 
over a one-kilometer path; the second method (recommended 
by the ACQUIRE’s User’s Guide) is to specify broadband 
atmospheric transmittance as a function of range. In the 
latter method the data may be obtained from measurements or 
predicted by using an atmospheric propagation model (e.g., 


LOWTRAN). In this thesis the SeaRad radiance model will be 


4‘] 


used for predicting atmospheric transmittance. Then outputs 
of this model will be included in the ACQUIRE data file. 

3. Sensor Model 

The ACQUIRE model permits sensor definition in two 
formats: either the sensor parameters are written in the 
ACQUIRE data file, or a separate sensor data file is 
included in the sensor look-up table. In both methods, 
required performance parameters (1.e., MRTD and MDTD), 
horizontal field-of-view (HFOV), and wide field-of-view - 
narrow field-of-view (WFOV-NFOV) ratio are included in the 
format. A data file name is required when using a look-up 
table in the main ACQUIRE file. This allows the user to 
define their own sensors and attach these user-defined 
sensors to the look-up table. 

4. Output Files 

ACQUIRE outputs are displayed in two formats. The 
first output file with the extension rl is automatically 
written for target discrimination performance after each 
tenet It lists the target discrimination ranges for 
different probability levels in a tabular form. The second 
one with the extension r2 is written when the user enters 
the appropriate command as explained in the User’s Guide. 


The format resembles the first one except that it lists the 
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a 


PanoermGn ser imMimcdt vom prebabilaties according to different 
detection ranges. The ri file is used to verify performance 
with respect to a specific probability requirement and 
“obtained by linearly interpolating the results listed in 
the r2 file, and, therefore subject to interpolation errors 
if the probabilities are changing rapidly with respect to 
the range increment.” [Ref. 13] 

ahiermregram Can alsemewrite the 12 OCuEDUE Eillew with “all 
headers and labels removed. Thus these files can be 
imported by plotting software and the outputs can be 


displayed in graphic form. 
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IV. COMPUTATIONS AND RESULTS 
A. SCENARIO INPUT PARAMETERS 

The scenario parameters were chosen to represent a 
beach environment for a joint operation. The parameters are 
related to air-to-ground weapon systems, which is the case 
that the TDAs under study deal with. These parameters were 
utilized to reach the ultimate goal of this thesis; that 
1s, to give ideas or make recommendations for a common TDA 
code that can be used in all services (Army, Navy and Air 
BoLce) . 

The target was chosen to be “Gunboat”, which is one of 
the targets in the WinEOTDA target look-up table. After 
studying its physical dimensions, it was noted that Gunboat 
has the same dimensions as R/V POINT SUR, which was used in 
the PREOS 92 Experiment in Monterey Bay in 1992. 

The atmospheric data were chosen to fit the properties 
of a typical operational environment of a naval target. An 
atmospheric data set collected from the Gulf of Oman, which 
was extracted from the EOTDA III model test data [Ref. 14], 
was used in this work. 

The published parameters of the SADA II Focal Plane 


Array (FPA) system were used to build a second-generation 


oa 


FLIR sensor using the FLIR92 model. The data relating to 
the physical parameters of the FPA were given by 
Ludwiszewski [Ref. 15]. The remaining input data required 
by the FLIR92 model were gathered from Ludwiszewski [{Ref. 
ie ad lp on second-generation sensor structure, and the 
textbook by Driggers et al [{Ref. 1]. The data set was used 
along with NETD and IFOV parameters as input to the FLIR92 
model to obtain MRT and MDT outputs. 

The following table summarizes the input data used in 


the WinEOTDA and ACQUIRE models. 


Target Gunboat 

Target Dimensions (m) Heo. 6 beer. Ss Wes ee 
Sensor User-defined SADAII 
Date/time (GMT) 07 Jan 199 syelg55 
Latitude 24deg 15min N 
Longitude 59deg 45min W 
Temperature (min) 

Temperature (max) 

Temperature (F) 

Dew Point Temperature (F) 

Sea Surface Temperature (F) 

Aerosol Model Index 6 (Navy Maritime) 
Vis toa der teams } S24 eda) 

Wind Direction (deg) 120 

Wind Speed (kts 

Low Cloud 


Inversion Height (ft) 
Table 4.1 - Scenario Input Parameters [Ref. 14]. 
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Be FLIR92 MODEL OUTPUTS 


The sensor model was formed using the input parameters 
included ina data file that was later saved in pure ASCII 
text mode. The model was run in the DOS environment for 
prediction of both MRTD and MDTD performance parameters, 
which would later form the input for ACQUIRE and WinEOTDA. 
THieGemeueopUrm Ores this @models comtaining the’ MRIDMand MbDTD 


predictions along with NETD and IFOV data is included in 


Appendix A. 


1. Comparison of FLIR92 Model Sensor Outputs With The 
Other Sensors in WinEOTDA Model 


The WinEOTDA sensors and the sensor built by the 
FLIR92 model were tabulated according to their physical and 
performance parameters. Later a comparison analysis between 
these sensors was performed to verify that the sensor built 
by FLIR92 had reasonable input and output data. The 
comparison charts are given in Appendix B, which includes 
it2ID, horizontal and vertical IFOVs and MRTD at min and max 
Spatial frequency comparisons. The sensor #127 represents 
the user-defined sensor, while the WinEOTDA standard sensor 
numbers range from 100 to 126. 

In the case of NETD it can be observed that the NETD 


parameter oof the user-defined sensor is within the 


oS 


theoretical limit, which is stated as 0.02 to 0.2 degree C 
by Driggers [Ref. 1], and in agreement with the other 
Sens omelets. 

Horizontal and vertical IFOVs, which are equal due to 
the square detector usage in the design process, are also 
reasonable when compared to the other sensor IFOVs. Since 
the user-defined sensor had a Narrow Field of View (NFOV), 
the comparison was made between the sensors of the WinEOTDA 
model. 

The maximum and minimum MRTD values showed an 
agreement between the sensors. However although the minimum 
Spatial frequency of the user defined sensor’s MRTD matched 
to the other sensor MRTDs, the maximum value was noticeably 
greater than those in the WinEOTDA model. This is because 
the sensor built by FLIR92 used the second-generation 
sensor (SADAII) parameters. Sensors of this generation 
offer increased resolution limits with a smaller detector 
size. Thus the MRTD values show an asymptotic behavior at 
higher spatial frequencies. 

The above comparisons allowed the conclusion that both 
physical and performance parameters of the sensor built by 
the FLIR92 model were reasonable and within theoretical 


ene ts. 
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ce SEARAD RADIANCE MODEL OUTPUTS 


The SeaRad Radiance model was used to calculate the 
atmospheric transmittance values required by ACQUIRE. The 
atmospheric data given in scenario input parameters were 
used as input to this model. The template and input data 
Given in Appendix C were used in the model, and 
transmissivities for O75) 0 ope tees Oo km ranges were 
obtained. An example output file along with a tabular and 
graphical form of all transmissSivities is shown in Appendix 
i) 

The maritime scenario used in this work represents the 
atmospheric conditions during winter in the Gulf of Oman, 
located in the sub-tropical region. Due to the location of 
the Gulf, even though the conditions were described for 
winter the mid-latitude summer model was chosen in SeaRad. 
The winter in the sub-tropics has almost the same 
atmospheric parameters as summertime in mid-latitude 


regions. 


D. WINEOTDA MODEL INPUTS 


The WinEOTDA model requires meteorological and 
tactical inputs before the model can be run. The 


meteorological data can be entered by selecting the MET 


2 


data summary From the maim menu™or @oyeeelickinicwecn the 
descriptor picture on the main screen. It includes the 
target location, surface weather characteristics at a 
specified time, and information about the boundary layer 
alienGg waszh cloud datax 

The tactical sLighalenaiv=yealiod Cie Operations and 
Intelligence (Ops) data include the inputs for _ sensor, 
target and backgrounds. The main menu or main screen can be 
used to enter the input data by selecting sensor, target, 
time over target or background. 

The detailed descriptions of the input parameters used 
in WinEOTDA can be found in Ref.12 and the WinEOTDA 
(Version 1.3.3) model’s help menu. 

1. Target Model 

Target and background information form this model. The 
target defines the size and physical characteristics used 
in WinEOTDA. Background data give information about general 
background, which is the dominant terrain feature of the 
target area, and the immediate area of the target. 

The target selected from the target menu ~was 
‘Gunboat’. For™ discrimimering the output ranges iE (@ Ne 
Gifferent viewing directions the heading, which refers to 


the direction of the target front aspect with respect to 
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north was entered as zero degrees. The operating state that 
gives the condition of the target at time over target (TOT) 
was selected ‘Off’, which meant that the target was heated 
by the environment alone. The entry form in Appendix E 
(Figure E.1) was used to enter the target data. 

The scenario conditions were described as a beach 
environment in the sub-tropical region. The background for 
this environment was selected as water, soil and vegetation 
respectively along with an ocean general background. Since 
the composition of the background has an effect on its 
heating and reflective properties, detailed information was 
entered for each background type by using the entry form in 
Appendix E (Figure E.2). 

The target location and time data are also used in the 
target model. They were entered according to the scenario 
input parameters. For time over target WinEOTDA offers two 
different options for decision makers and operators. The 
model can calculate the output ranges for either the 
execution phase or the planning phase of the scenario. 
While the execution phase needs the exact operation time 
and computes the output ranges according to varying viewing 
directions, the planning phase uses time intervals to 


calculate the detection ranges. Both parameters can be 


oe) 


entered by using the TOT entry form in Appendix E (Figure 
E.3). The execution phase was chosen in this thesis for 
comparison purposes of WinEOTDA to ACQUIRE. 


2. Atmospheric Model 


The atmospheric model uses the input Meteorological 
and Site (Met) data to calculate the atmospheric 
transmittance. The Met information is entered using the 
input form in Appendix E (Figure E.5) via the main menu or 
the main screen. The entry form offers more options that 
can be seen after clicking on the individual parameter 
labels. Furthermore, the graphical view of each entry can 
be displayed by ranlie fare, oll C kaiacia ore the indasgaeiaa i: 
parameters. The weather forecast data can be entered as a 
spot entry, which is the case in this thesis, or as 24-hour 
cycle data. 

The following set of meteorological data from the 


scenario parameters was used in WinEOTDA model. 
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O7evian 1993, 1055 
24deg 15min N 


59deg 45min W 









(GMT) 





Date/time 
Latitude 
Longitude 
Temperature (min) 
Temperature (max) 
Temperature (F) 
Dew Point Temperature (F) 
Sea Surface Temperature (F) 
Aerosol Model Index 

Were ophlabiey am [Oily 

Weather 

Wind Direction (deg) 

Wind Speed (kts) 
Low Cloud 



























6 (Navy Maritime) 
V5 ei 2A ai) 










120 














Height 
Amount 













Inversion Height (ft) 
Upper Level (UL) 

UL Temperature (F) 
UL Dew Point Temperature (F) 
UL Aerosol 

Wie visibility 
Table 4.2 - Meteorological Input Parameters [Ref. 14]. 









3. Sensor Model 

The sensor entry form that can be reached via the main 
screen or the main menu, included in Appendix E (Figure 
E.4), offers a numbered sensor list to the operator. The 
user-defined SADATI second generation FLIR sensor was 
attached to this list as sensor #127. The NETD, XIFOV and 
PEEey, MRTD and MDTD Gutputs of the FLIR92Z model formed the 
sensor data file. 

The remaining parameters for the entry form include 


the sensor height, viewing direction and scene complexity. 
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The sensor height was entered as 500, 2000, and 4000 ft 
respectively, along with a zero degree viewing direction. 
The scene complexity input describing the number of objects 
in the immediate vicinity that can be mistaken for the 


targets was chosen to be ‘None’. 


EB. WINEOTDA MODEL OUTPUTS 


The scenarios as previously defined were run for each 
of three different sensor altitudes: 500, 2000, and 4000 
fires After each “successful. runiesateandicaticon sci aa 
‘successful run’ was presented at the bottom of the screen 
and the maximum detection ranges for different targets with 
different backgrounds were displayed in tabular form on the 
main screen. The other output files eRe od a Chapi-ei 
III were also automatically generated. The examples of 
these output files along with a main screen output table 
are given in Appendix F. The alphanumeric output summarizes 
some of the input parameters and displays the calculated 
Detection Ranges, Thermal Contrast, and Target Temperature 
for the permutation of targets #1 and #2 and three input 
backgrounds, where target#l and target#2 are the same. The 
same results are also displayed in graphical (Figure F.1) 


and tabular (Figure F.2) formats along with a main screen 
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output (Figure F.3). These outputs were obtained with the 
original scenario input parameters for the sensor at 2000 
ft altitude. 

In these scenarios the heading of the target and the 
viewing direction of the sensor were chosen to be zero 
degrees. This information was utilized to determine the 
front and side aspect angles in the output files as: 000 
for Front and 090 for Side views. Although the model gives 
outputs for 45 degree intervals in viewing direction, only 
the 000 and 090 degree directions were used to compare the 


results with the ACQUIRE outputs. 


F. ACQUIRE MODEL INPUTS 


The ACQUIRE model uses a data file which includes 
target, sensor, and atmospheric information to predict the 
discrimination ranges. This file must be in pure ASCII code 
and written before the model is run. The data file format 
can be seen in Appendix G, where the section between the 
header line and ‘end’ of the file forms the input data file 
format. 

1. Target Model 

The ACQUIRE has two different ways of defining a 


target in the model: The target look-up table, which 


oil 


consists of a number of land targets, can be used, or user- 
defined data can be entered. The latter was used to specify 
the size of the naval target, Gunboat, which was not 
represented in the internal look-up table of the model, and 


the following format was used to enter target data: 


Lace ec & 
characteristic _dimension 0.0 meters 
Signature 0.0 degrees_C 


The signature parameter, which 1s the temperature 
difference between the target and its surrounding 
background, was taken from Shumaker et al [Ref. 10] where 
the representative values of ship differential temperatures 
were readily available for summer conditions at the 
specified operation time eo Os ie AM). The selection 
procedure from the same reference was followed and the 
Signature was determined to be 7.71 degrees C. 

The selection of a proper target size, represented by 
“Characteristic cimension - ts escritical. For discriminat tom 
the target size directly affects the number of spatial 
cycles that can be resolved across the target. In the 
ACQUIRE User Guide a characteristic dimension that is equal 
to the square root of the projected area of the target in 


Mewees 1S recommencdeamior caleculations. ms mecessi_atces 
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the proper calculation of the projected area of the target. 
The program allows alternative definitions for target- 
projected area to be applied Ee off-menu targets. For this 
thesis the following equations given by Shumaker [Ref. 10] 


were used in projected area (A,) and critical dimension 


(D.) calculations with the following target orientation 


c 


model: 


LOS | 
NN 





Figure 4.1 - Orientation of Targets “After [Ref. 10]”. 


A, =lhcos@cos¢@+whcos@sin @+lwsin@ (ay) 
where 
l 1s the actual target length (m) 
W 1s the actual target width (m) 
h 1s the actual target height (m) 
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g 1s the elevation angle (deg) 


Q 1s the azimuth angle (deg) 


D, = {Ar (Aes 


Due to the importance of the subject the following 
section will cover an analysis of critical dimension 
caleulaeron . 

a) Critical Dimension Analysis 

The ACQUIRE code uses a critical dimension 
parameter calculated off-line for Nee eeee ied targets. 
This value is entered as a constant in the model and is not 
calculated continuously within different time intervals for 
varying ranges and altitudes. 

The orientation model in Figure 4.1 and the Equations 
4.1 and 4.2 were used to calculate critical dimension 
values for different ranges (between one and 30 km) and 
altitudes (500, 2000, and 4000 ft). The altitudes were 
chosen to represent a descending aircraft carrying the 
sensor. 

The calculations in Appendix H showed that the 
critical dimension is very similar for longer ranges at 
different altitudes, but at short ranges the results for 


different altitudes are dramatically different. Varying 
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ranges and altitudes definitely affects the aspect angle, 
which causes a change in apparent target size. The critical 
dimensions were calculated for front (F) and side (S) views 
of the targets, which are represented by the azimuth aspect 
angle ‘90’ and ‘OQ’ degrees. 

To evaluate the impact of the changes in aspect ratio 
due to altitude changes, critical dimension values for 500, 
2000, and 4000 ft sensor heights and 25 km range were 
calculated. These were used as input parameters for the 
target model. 

2. Atmospheric Model 

The SeaRad Radiance model outputs in Appendix D were 
used in building the ACQUIRE data file. This is the 
recommended method in the ACQUIRE User’s Guide gene 
calculating atmospheric transmittance and the following 
format is the only way of entering these data: 

>band_averaged_atmosphere 


# points: 0O km transmittance 
The km array holds the range in kilometers and 


requires at least three and not more than fifty points with 


the data ordered by increasing range. The transmissivity 
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array holds atmospheric transmittance with the same number 
of points as the km array. 


3. Sensor Model 


ACQUIRE models the sensor in two ways as mentioned 
earlier in Chapter III. In this work the second method was 
adopted in which a separate sensor data file that is later 
included in a sensor look-up table was used. The outputs of 
the FLIR92 model, which include horizontal, vertical, and 
2D MRTD values, were utilized to build the sensor data file 
with the following format: 

Header line 

>systemA 

WEov: O20 viov: OF0SwWwmtew ratio: 0.0 
@MRTD_2d 


HemOIntS. CU eR g phar. (ra eeemebe MRTD 


The first line of the file is a ‘Header line’ that is 
always required. The second line is the name of the sensor 
system and remaining lines contain sensor physical and 
performance parameters. For discrimination purposes 2D MRTD 
outputs of the FLIR92 model were used in this thesis. This 


externally built sensor data file was included in the 
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ACQUIRE data file to determine the sensor system with the 


fel 1 OwamGmer ote matey 


>sensor_lookup 


data_file_name : 1s the name of the file in 
which sensor performance curve data is stored in 
lookup tables. 

sensor_id ac Slices = ldkbela — sValsienqu tee 
select from the lookup table. 

performance_mode: selects the performance data 


to read from the lookup table. 


The target discrimination criteria listed in Table 2.2 


in Chapter II were entered using the following format for 


the MRT ‘performance_mode’ of the sensor model: 


>cycle criteria 


WFOV 


qetection 


detection _n50 0. 0s ferme WEev 
detection_n50 O20 Sapo oN OV 
classification _n50 C2 Oe tor troy 
recognition _n50 OG te TONE Oy 
PASH ELTICAation 150 CSO: for NERO 

is only used for detection purposes. After 


of the target NFOV is used to resolve the details 


nee AtSserininat lon DUnDOeseS - 
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G. ACQUIRE MODEL OUTPUTS 


The input parameters including target, atmosphere and 
sensor data along with the varying characteristic size of 
the target according to different altitudes were entered 
into the following ACQUIRE data files: 

e SadallIfa: is the ACQUIRE data file for front view 
of Gunboat seen from the sensor at 500 ft altitude 
and 25 km range, 

e Sadallifb: 1s the ACQUIRE data file for front view 
of Gunboat seen from the sensor at 2000 ft altitude 
and 25 km range, 

e Sadalilfc: is the ACQUIRE data file for front view 
of Gunboat seen from the sensor at 4000 ft altitude 
and 25 km range, 

e SadalIsa: is the ACQUIRE data file for side view 
of etnboat seen from the sensor at 500 ft altitude 
and 25 km range, 

e SadalIsb: is the ACQUIRE data file for side view 
of Gunboat seen from the sensor at 2000 ft altitude 


and 25 km range, 
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e SadalIsc: Tomei “ACOUlnM "Gata Lile for side vlew 
of Gunboat seen from the sensor at 4000 ft altitude 
and 25 km range. 

The output files generated after each run of the above 
PIles= with  themwextensiton “Of “rl and =12 “ase imcluded in 
Appendix I. Since the model uses the same file name for the 
output as the input, all output files have the same names 
as given above. The r2 extension files provided input for 
Windows excel toe produce =the plots ofmthe” probaomirty or 
discrimination given range parameters. These plots were 
also included in Appendix I along with the alphanumeric 
results. 

As can be seen from the output file formats, they 
basically use three different summary sections to display 
mime results. The first is the one containing intemmatvon 
about the input ACQUIRE data file used, the second is the 
messages containing information about the sensor structure 
and a detailed examination of intermediate results. And the 
last section gives the discrimination ranges versus’ the 


Given probabilities. 
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H. COMPARISON OF WinEOTDA AND ACQUIRE MODEL OUTPUTS 


The TDAs under study were compared under three 
arveiterens conditions. Pade Sits the Sfrginal scenario 
parameters previously listed were used for range prediction 
in both codes. Second, WinEOTDA Met data were modified to 
get the same 4km transmissivity value as in ACQUIRE. Last, 
the two codes were compared fe ©) 1G the same 4km 
transmissivities obtained by using the Beer’s law 
approximation in the ACQUIRE model. This section will cover 
the procedure and the results related to the first 
Sonat von: 

The models were run for different sensor altitudes 
(S00, 2Z000h) amc 400GEt) with Chemortcqinal sscenariomanpuE 
parameters. The following results (Table Lo) were 


obtained. 


WinEOTDA DETECTION RANGES ACQUIRE DETECTION RANGES 
(4 KM TRANSMISSIVITY = (4 KM TRANSMISSIVITY = 0.407) 


SENSOR 0.60) 


— SIDE VIEW FRONT VIEW SIDE VIEW FRONT VIEW 
= deg) (0 deg) (90 deg) (0 deg) 


Table 4.3 - WinEOTDA and ACQUIRE Detection Range Comparison 
Table With The Original Scenario Parameters for Different 
Sensor Altitudes 
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The obvious difference in Table 4.3 was found to be in 
4km transmissivities calculated by two programs. This 
difference is due to the detailed input data structure of 
WinEOTDA model as opposed to ACQUIRE. Two crucial 
parameters, temperature and dew point temperature, io 
calculating relative humidity, which 1s an input (for 
atmospheric extinction coefficient calculations, are not 
entered by the users in ACQUIRE atmospheric model (i.e. 
SeaRad in this study). Instead the atmospheric 
transmittance model uses the default values for the 
specified region. Additionally, the low cloud cover, which 
was ayo used in ACQUIRE, caused a difference sal 
predictions. 

The results in Table 4.3 showed that the two programs 
give different detection ranges with the same scenario 
mou t parameters. WinEOTDA predicted a longer-range 
performance ee a1 for NFOV for both aspects. It output 
almost twice as long detection range as ACQUIRE at 500ft 
sensor altitude in NFOV detection. This was not observed 
for WFOV detection. 

As seen in Table 4.3 WinEOTDA detection ranges vary 
with the changing sensor altitudes, but this is not 


observed for the ACQUIRE model, or at most the observed 


qa 


change 1s very small. This is because of the difference in 
target modeling in the two programs. While WinEOTDA uses a 
powerful target signature model, TCM2, which calculates the 
mean temperature and projected areas within specified time 
intervals for a user defined period, ACQUIRE requires only 
a measured target signature and critical dimension value 
for off-menu targets, ignoring the changing sensor altitude 
effects. However, as analyzed before, the varying sensor 
height affects the projected area of the target and 
consequently critical dimension. In WinEOTDA, calculated 
total projected area is used in determining the sensor 
performance model detection ranges as a function of target 
Spatial frequency. 

On the other hand, background information which is 
useful for target signature calculation is ignored for user 
defined targets in ACQUIRE. Despite the convenience of 
describing off-menu targets in ACQUIRE by their measured 
parameters, this method depends heavily on the 
accountability of the reference used. However WinEOTDA lets 
the users enter a detailed background input data set that 
is used by the TCM2 model for target signature calculation. 
Thus the model can generate a more accurate target 


Signature. 
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Es COMPARISON OF WinEOTDA AND ACQUIRE OUTPUTS WITH THE 


MODIFIED WinEOTDA METEOROLOGICAL DATA 


The previous section showed that the models calculated 
different 4km transmissivity values for the same original 
scenario input parameters. In order to see whether the same 
transmissivity values will give the same results, some 
modifications were made in the WinEOTDA meteorological 
input data. Temperature and dew jsrevsauc temperature 
parameters were changed until the same 4km transmissivity 
as in ACQUIRE was found. After achieving the same value, 
Ene | program wase sun “Eor 5007, 2000, and 4000 SEe sensor 
altitudes and the results in the following table were 


@otained. 











ACQUIRE DETECTION RANGES 
(4 KM TRANSMISSIVITY = 0.407) 


WinEOTDA DETECTION RANGES 
(4 KM TRANSMISSIVITY = 
Oran} 

(90 deg) (0 deg) (90 deg) (O deg) 
Table 4.4 - WinEOTDA and ACQUIRE Detection Range Comparison 
With The Modified WinEOTDA Meteorological Parameters To Get 
The Same 4km Transmissivity Value As In the ACQUIRE for 

Different Sensor Altitudes 















SENSOR 
mele Tl 
at ) 







Table Ae shows that even qoge the same 4km 


transmissivities are used in both models, the predicted 
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detection ranges are still different from each other and 
WinEOTDA predicts longer-range performance For NFOV 
detection for both aspects. Since WinEOTDA still uses the 
Beer’s law approximation in calculating transmissivities 
for the other ranges, this 4km transmissivity will be the 
only common value Gog WinEOTDA and ACQUIRE. SUguLS 
approximation along with a better target signature and 
critical dimension calculation in WinEOTDA cause _ the 
difference in the output ranges. 

As seen from Table 4.4, the WinEOTDA detection ranges 
are reduced, as expected, when compared to those in Table 
4.3. This 1s because the atmospheric model uses a lower 4km 
transmissivity as a reference to calculate the 
transmissivities ies different | ranges. The smaller 


transmissivities cause a decrease in the predicted ranges. 


J. COMPARISON OF WinEOTDA AND ACQUIRE OUTPUTS WITH THE 
USE OF BEER’S LAW APPROXIMATION IN ACQUIRE MODEL 


The predicted detection ranges of the WinEOTDA and 
ACOULRE models were compared iets the same 4km 
transmissivities as in the previous section. But this time 
the Beer’s Law approximation was also used in the ACQUIRE 


atmospheric model. First the following table was built 
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using the ACQUIRE detection ranges in Table 4.3 and the 
outputs obtained after running the ACQUIRE model for 
different sensor altitudes (500, 2000, and 4000ft) with the 
transmittivities calculated by the Beer’s Law 


approximation. 












ACQUIRE DETECTION RANGES in 

Km (BEER’S LAW APPROXIMATION 

FOR 4 KM TRANSMISSIVITY = 
0.407) 


(90 deg) (0 deg) (90 deg) (0 deg) 


Table 4.5 - ACQUIRE Detection Range Comparison With The 
SeaRad Atmospheric Transmittance Parameters and Modified 
ACQUIRE Atmospheric Transmittance Parameters by Beer’s Law 
Approximation for Different Sensor Altitudes. 


ACQUIRE DETECTION RANGES in 
Km (4 KM TRANSMISSIVITY = 
0.407) 















SENSOR 
HEIGHT 
it) 








Then the predicted ranges in Table 4.5 obtained by 
using Beer’s law were utilized to compare the results of 
the two models for the same atmospheric model inputs. Table 


4.6 displays the results of this procedure. 
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WinEOTDA DETECTION RANGES ACQUIRE DETECTION RANGES 
(4 KM TRANSMISSIVITY = (BEER’S LAW APPROXIMATION FOR 


Sra 0.41) 4 KM TRANSMISSIVITY = 0.407) 


ec cie SIDE VIEW FRONT VIEW SIDE VIEW FRONT VIEW 
(90 deg) (O deg) (90 deg) (0 deg 


| 500 | 31.9 | 18.4 | 29.9 


Table 1.6 - WWinEOTDA ra ACQUIRE Det earn Range Conpari aon 
Table With The Modified ACQUIRE Atmospheric Transmittance 
Parameters To Get The same 4km Transmissivity Value As In 
WinEOTDA (Beer’s Law Approximation) For Different Sensor 

Altitudes 





As seen from Table 4.5, the Beer’s law approximation 
reduced the range performance of ACQUIRE. This 1S an 
expected result, since the exponential Beer’s law gives 
smaller transmissivity values than the SeaRad model. A 
MathCAD file was included in Appendix J to show that the 


Beer’s law approximation gives smaller transmissivities 


than SeaRad model where 7) represents the transmissivities 


calculated by using the Beer’s law approximation and To,oraq 


displays the outputs of SeaRad Radiance model for different 
ranges. 

Table 4.6 shows that WinEOTDA predicts better 
performance for both aspects in NFOV and WFOV with the same 
atmospheric model inputs. The difference is due to the 


better structure of the WinEOTDA target model to calculate 


m6 


the target signature and better approximation fe © 12 
Calculating "CrittealwcemensioOn of the target. These two 
parameters were frequently mentioned in the last sections 
of this thesis as the main causes of differences in 
predictions. However WFOV detection ranges do not vary 


according to the changing transmissivity values. 
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Was CONCLUSIONS AND RECOMMENDATIONS 

This thesis has presented an analytical comparison 
between the Army FLIR TDA, ACQUIRE Version 1 dated 1995, 
and the infrared module of the Navy/Air Force TDA, WinEOTDA 
Version 1.3.3 dated 1998. The programs were compared with 
respect to different means they used to model target, 
atmosphere and sensor. They were analyzed for the same 
scenario conditions, in which the scenario parameters were 
chosen to reflect a beach environment within the GOnecepe OL 
‘Joint operations’. 

The research questions addressed were to find the 
differences in the modeling of gounderlying physical 
peinciples, input parameters, and predicted detection 
ranges; suggestions were sought for modification of the 
codes that would lead to equivalent output for the same 
inputs. The following sections will give responses to these 
questions that may finally determine the possibility of 
uSing one of the codes under study as a standard TDA for 


ome SOY Vices. 
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Po COMMON SET OF INPUT PARAMETERS 


The common set of input parameters that can be used to 
operate both WinEOTDA and ACQUIRE can be grouped as target, 
atmosphere and sensor data. As studied in Chapter III these 
data form the inputs for the models, which are known by the 
same names, constituting both TDA codes. Although both 
programs have the same model structure, the treatment of 
inputs shows some differences that will be made clearer in 
the following sections. 

The target models in this work were built using a 
naval target and its related backgrounds. The R/V POINT 
SUR, which happened to have the same dimensions as the 
Gunboat entry in WinEOTDA target look-up table and 
backgrounds depicting a beach environment formed the inputs 
for both models. While no external calculations for target 
Signature and critical dimension were needed in WinEOTDA, 
ACQUIRE required these computations for this off-menu 
Earger 2 eeUtiimes internaloommp table 1s wsetuletor aime 
targets only, while WinEOTDA can be directed at maritime 
and overland scenarios. 

The atmospheric data were entered according to the 


scenario parameters. WinEOTDA required specific parameters 
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for its internal transmittance calculations. On the other 
hand in ACQUIRE these computations were externally handled 
by the user and the results were incorporated into the 
program. 

The sensor models in both codes utilized the FLIR92 


model outputs. 


B. DIFFERENCES IN THE PREDICTIONS OF CODES WITH THE 

COMMON INPUT PARAMETERS 

Using equivalent data the two programs yielded 
different detection ranges for different sensor altitudes. 
Table 4.3 in the previous Chapter provided these detection 
ranges of both codes. While varying sensor altitudes caused 
only a slight change in ACQUIRE detection ranges, WinEOTDA 
displayed significant differences for the same sensor 
altitudes. 

WinEOTDA predicted a longer detection range than 
ACQUIRE 100% of the time. In particular, NFOV detection 
ranges of WinEOTDA were twice as long as the ranges in 
ACQUIRE. However, although the WFOV detection ranges of 
WinEOTDA displayed better performance than ACQUIRE, they 
seemed to be insensitive to differing aspect angles and 


showed a very small change for different sensor altitudes. 


on 


For ACQUIRE, WFOV detection ranges did not change as much 
with changing sensor altitude as did the NFOV detection 


ranges. 


Cc DIFFERENCES IN THE TREATMENT OF CODE INPUTS 


WinEOTDA and ACQUIRE have the same input models. 
However, there are some differences in the treatment of 
inputs, and these differences in the interaction of sensor 
and target give rise to the differences in predictions. 

The WinEOTDA target model uses TCM2 to calculate the 
target signature. This model treats the target as a mesh of 
different nodes and calculates the resultant signature by 
considering the transfer of heat between these nodes. It 
also takes into account the atmospheric effects on the body 
heat of the target surface and calculates the signature 
according to the input meteorological data. It is believed 
to be the most accurate model available to calculate the 
target signature, and the results are certainly better than 
the ACQUIRE target model. The calculated target-to- 
background temperature difference at zero range will have a 
great effect on the discrimination of the targets. On the 
other hand, the ACQUIRE target model allows the user to 


input the critical dimension of the target and the measured 
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target-to-background temperature difference as the ‘target 
Signature’ for off-menu targets. 

There are some deficiencies in this method. In the 
case of target signature, the difference between the 
target’s average surface temperature and the average 
immediate background temperatures at zero range are used 
for calculations and it cannot be as accurate as the TCM2 
Model used in WinEOTDA. Further, only one value of critical 
dimension is required. However after the analysis of the 
critical dimension calculation by using aspect angles, it 
was noted that even though the impact of the change in 
altitude for longer ranges has a negligible effect on the 
critical dimension, it has a considerable effect at shorter 
ranges and the critical dimension changes significantly for 
shorter ranges at different altitudes. Although the longer 
ranges were used for critical dimension calculation in this 
thesis, varying sensor altitudes lead to differences in 
range predictions. As can be seen in Appendix H_ the 
Sritical sdimensmencmatweeSeekm for sa00, 000, amd 4000 ft 
altitudes are very close but not identical. These small 
differences cause the changes in detection ranges for 


varying sensor altitudes. 
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Due to the unavailability for analysis of the original 
ACQUIRE source code, the treatment for the targets included 
in the look-up table could not be examined. Thus comparison 
of the treatment of target models in the two codes could 
only be made for external targets in ACQUIRE. The author’s 
suggestions are based on these results. 

The WinEOTDA atmospheric model evaluates atmospheric 
transmittance values using a limited version of the LOWTRAN 
model according to the meteorological data input by the 
user, which implies that more detailed and accurate data 
can give better outputs. The transmittance value at 4km is 
calculated by LOWTRAN and then the extinction coefficient 
obtained for this range is used to calculate the other 
transmittance values using the Beer’s Law approximation. 
But the use of the broadband Beer’s law approximation is 
known to give erroneous results as discussed in Chapter II. 
Thus, this must be considered as a weakness of this code. 

The range predictions of the two programs were 
compared for the same 4km transmittances obtained by 
modifying the meteorological data of WinEOTDA. But the 
Beer’s law approximation was still used in WinEOTDA as 
opposed to SeaRad transmittance used in ACQUIRE. It was 


observed that the WinEOTDA NFOV detection ranges decreased 
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and got closer to ACQUIRE predictions. This change was due 
to the modification ube! temperature and dew point 
temperature parameters of the WinEOTDA meteorological data. 
Temperature was reduced and dew point temperature was 
increased, which caused an increase in RH. But the WFOV 
detections were unchanged and were insensitive to changing 
altitudes and transmissivity. Amother significant result 
was the diminishing detection range with increasing sensor 
altitude in WinEOTDA NFOV detection. Here the only 
meteorological parameter modified was dew point temperature 
to get the desired change in relative humidity to match 
transmissivity between the models. This is thought to be 
the only effect that caused the change. 

The ACQUIRE atmospheric model allows two different 
methods for calculating the atmospheric transmittance 
values. In the first method, the transmittances are 
calculated by the broadband Beer’s law approximation, which 
has the same deficiencies as mentioned for the WinEOTDA 
model. In the second method, recommended by the ACQUIRE 
User Manual, one of the atmospheric transmittance codes is 
used to get the transmissivities. Then the resultant 
transmittance values are directly input into the ACQUIRE 


fiaetertlow Tae oademor thinks that this treatment is better 
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than the first and the method used in WinEOTDA. The second 
method was used in this thesis to get the required 
atmospheric transmittance values. Furthermore, different 
effects of the first and second methods on the model 
outputs were analyzed as described in Chapter IV. 

ACQUIRE model atmospheric transmittances were 
calculated by using the Beer’s law approximation as in 
WinEOTDA. For the same 4km transmissivity values and method 
for calculating the transmissivities for the other ranges,. 
the detection ranges of both programs were compared. The 
outputs are shown in Table 4.5 and Table 4.6 in Chapter IV. 
The Beer’s law approximation used in ACQUIRE reduced the 
detection ranges as expected. This caused a _ greater 
difference between the two codes’ NFOV detection ranges. In 
the case of WFOV detections the insensitivity of WinEOTDA 
to varying sensor altitudes and aspect angles was still 
observed. 

WinEOTDA sensor model utilized the FLIR92 outputs to 
build a sensor data file. The same procedure was also 
followed by ACQUIRE. Although the input parameters and 
sensor data file structure are similar in both codes, 


ACQUIRE predicts classification, recognition and 
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identification detection ranges in addition to WinEOTDA’s 


particular detection range predictions. 


D. SUGGESTIONS FOR THE MODIFICATION OF CODES 


Some suggestions can be made about the modification of 
the codes that would lead to equivalent output for the same 
inputs. Firstly, since the FLIR92 model was commonly used 
to create a sensor data file for both codes and the sensor 
models have the same structure except for their different 
ways of handling the sensor data, FLIRG2Z can be 
incorporated into both programs for modeling the new 
sensors. In fact ACQUIRE currently uses the outputs of the 
FLIR92 Model, but WinEOTDA can also be integrated with 
FLIR92 for building new sensors either externally or 
internally. However the cost and time must be taken into 
consideration before the integration. 

The WinEOTDA target model uses TCM2, which is seen as 
the best and the most accurate model available for target 
Signature calculations. Thus, the only suggestion for the 
modification of target model might be to include more 
targets in its look-up table. On the other hand, since the 
ACQUIRE program code was not available for examination, the 


treatment of its target model for look-up table targets is 
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unknown and the suggestion will be to include more targets 
in the model. However for off-menu targets, a target 
Signature model such as TCM2 can be included to get better 
and more accurate predictions. 

Atmospheric transmittance values required by the 
WinEOTDA atmospheric model can be calculated completely by 
LOWTRAN/MODTRAN and used to find the apparent temperature 
difference of the target-to-background. This is expected to 
give more accurate results in the prediction of the model. 
However a trade-off analysis must be performed before the 
integration of the whole model, as in the sensor modeling 
case previously mentioned. The cost and time needed to run 
the program will increase and this will cause some 
problems. Especially when the importance of minimizing the 
time required to reach a decision is considered, it will 
not be easy just to decide on the integration of the whole 
program before an exhaustive analysis. 

In the ACQUIRE atmospheric model use of the second 
method is recommended to avoid the erroneous results of the 


broadband Beer’s law approximation in the predictions. 
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E. CODE SELECTION FOR INTER-SERVICE USE 


The WinEOTDA code seems to have a deficiency in its 
atmospheric model, which is not easy to fix for the reasons 
given in the previous section. But the target and sensor 
models are powerful and give accurate results. 

On the other hand, ACQUIRE has some shortcomings in 
modeling targets and backgrounds, and in the method for 
transmittance calculation. Furthermore it 1s not user 
friendly, and requires some codes to be written in specific 
formats to run. It also requires an operator trained in IR 
mNecOry and@mehe “operaeion Mot the code, Swhiche is Snot 
generally available na! the operational environment 
envisaged for naval TDA use (i.e. ordnance selection and 
pre-sortie mission planning.). 

Although the range predictions of the two programs 
compared in this work have not been validated by real world 
measurements, the better performance of WinEOTDA, its easy 
to use structure and powerful target model display an 
advantage in choosing a standard TDA for inter-service use. 
However the author believes that at the unclassified level 
without using the real sensor data and predicted ranges, it 


1s not easy to decide on a standard code. 
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F. RECOMMENDATIONS FOR FURTHER RESEARCH 


The differences in the predicted ranges as shown in 
Chapter IV point out the need for field testing of the two 
programs to determine the accuracy of detection ranges. If 
the same results occur ina field test, one of the programs 
may show as better than the other. This could result in an 
improvement to the other program, or choosing the better 
one as a standard TDA for inter-service use. 

The comparisons on an unclassified level might not 
reflect the actual performance of the Ses The real 
sensor parameters and predicted detection ranges can be 
more useful to prove the reliability of performance. Thus a 
classified level research study with all the needed real 
world parameters will give better information to decide on 
or mModtTiy@a specific TDA code. This would require a 
measurement campaign on the level of the MAPTIP [Ref. 20] 
or EOPACE [Ref. 21] international measurement series. 

The next level of comparisons must take place between 
ACQUIRE and the Target Acquisition Weather Software (TAWS), 
which 1s an upgrade to the EOTDA program. The summary 
information about TAWS and the comparison tables of delta T 


and detection ranges of WinEOTDA and TAWS can be found in 
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Appendix K. A more systematic and detailed comparison of 


the two codes is recommended for future study. 
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APPENDIX A. FLIR92 MODEL OUTPUTS 


U.S. Army CECOM NVESD FLIR92 The Meweeewe 5 24:55 2000 
eutput file: SADAIZI.1 Sigm ears cing 

data file: SADAII 

command line arguments: -d SADAII -o SADAII -p MRT 

begun data tile @aisting ~ .... 

gen2: sample data file for 2nd generation FLIR with SADA II FPA 


>spectral 
spectral_cut_on 30 microns 
spectral_cut_off a microns 
adiffraction_wavelength 020 microns 
SOODLICS  - 
f_ number G0 -- 
eff focal_length ZOzAO em 
eff_aperture_diameter 10.0 cm 
optics_blur_spot C201 mrad 
average_optical_trans Cray -- 
SODtLIcCS—2 
HFOV:VFOV_aspect_ratio C.0 -- 
magnification 0.0 -- 
frame_rate 2020 Hz 
fields_per_frame le @ -~ 
>detector 
horz_dimension_(active) 25.24 microns 
vert_dimension_(active) 2a oA microns 
peak_D_star 1.5e10 em-sqrt (Hz) /w 
integration_time 02007 microseconds 
1/f_knee_frequency 30 Be 
>fpa_parallel 
# detectors_in_TDI a) -- 
#_vert_detectors 480.0 —- 
#_samples_per_HIFOV Ze -- 
#_samples_per_VIFOV 250 -- 
3d0B_response_frequency 0,5 Pale, Hz 
scan_efficiency Oi es -- 
>electronics 
high_pass_3db_cuton ILE) Hz 
high _pass_filter_order eLaie) -- 
low_pass_3db_cutoff LOCOCO Cs.6 Hz 
low_pass_filter_order C20 -- 
boost_amplitude 0.0 -- 
boost_frequency 0.0 Hz 
sample_and_hold HORZ NO_HORZ_VERT 
>display 
display_brightness 1LIO Se, milli-Lamberts 
display_height 15.22 em 
display_viewing_distance 50.0 cm 
>crt_display 
#_active_lines_on_CRT 480.0 -- 
hOuzZ2crt spot_sigma OF 16 mrad 
vert_crt_spot_sigma 0.0 mrad 


oS 


>eye 


threshold_SNR 25 -- 
eye_integration_time opr sec 
MTF EXP Bor _or_NL 
>random_image_motion 
horz_rms_motion_amplitude O02 mrad 
vert_rms_motion_amplitude 0202 mrad 
>Sinusoidal_image_motion 
horz_rms_motion_amplitude ome 0, mrad 
vert_rms_motion_amplitude 0.0 mrad 
>3d_noise_default 
noise_level MOD NO_LO_MOD_or_HI 
>spectral_detectivity 
# points: 9 microns detectivity 
8.00 02666 
8.50 0.708 
9 2010 Omr'5 0 
ole OS 2 
O20 0 0.833 
10.50 0.875 
i ie O60 Geo17 
ull) 0, 0.958 
12 00 iO 
>end 
end data file listing 
MESSAGES 
dtagnostic(): Using default@3bp noise compoenentse 
diagnostic(): Using _MOD_ level 3D noise defaults. 
diagnostic(): Diff. wavelength set to spectral band midpoint. 
dtaqostuei): HFOY:-VFOV aspect ratio defaultedera 1.33. 
diagnostic(): Fields-of-view calculated by model. 
dagmeseic(): Electronics high pass filter Gdetaulted to cormdexrgl. 
diagnostic(): Electronics low pass filter defaulted to order 1. 


CALCULATED SYSTEM PARAMETERS 
field-of-view: 2.323 oo 


20. SSRs O24 Sy 

Magniticarien: 623.23 
Optics bling coe: 48.800 
0.244 

detector IFOV: O.127n eer, 
scan veloctuy. 62129 
dwell time: 7.833e-005 


TEMPERATURE DEPENDENCE 


parameter NETD @ 300 K 

white NETD Of135 deg C 0.000 
classical NETD Ole 5.deq7e€.0.000 
sigma_TVH NETD O-Os deg =c 0.000 
Sigma_TV NETD O2077 sdeg e.0. 006 
Sigma_V NETD O07 7 edegee..0- 000 


degrees 

mrad 

microns (diffraction—-mmmred) 
mrad 

mrad 

mrad/second 

seconds 

NETD @ OK noise bandwidth 
deg C 1.003e+004 Hz 

deg C 1.007e+004 Hz 

deg C 3.134e+003 Hz 

deg C 

deg C 
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Planck integral 1.978e-004 0.000e+000 W/ (cm*cm*K) 


w/D-star 2.439e+006 CeOQO00e- 000) seme (Hz)/(cm*K) 
PREP ISIER VALUES AT NYOUTS?. 
horz BePRE (7.07) =s02000 vert H_PRE(7.87) = 0.000 


SAMPLING RATES 
horizontal 15.75 samples/mr 
vertical 15.75 samples/mr 
effective 15.75 samples/mr 
SENSOR LIMITING FREQUENCIES 


spatial Nyquist 
horizontal yes 7 Toe84 
vertical Tao TaSe 
effective Paod 7.87 

MRTD 3D NOISE CORRECTION (AVERAGE) 

300° K OK 
Horizemzal, 2.000 G0 00 
vertical 32263 O77 000 


MRTD AT 300 K BACKGROUND TEMPERATURE 


CV/mr Nee ey /mz vert ey / Tas 2D 
05 0.394 0.007 C. 0s 0.394 0.065 0.830 0 
mo. 060.787 O20 ne Oe i(0, C2787 OF 100 ge 28 0 
telo) 1.181 0.0352 0.55 aged spa we. 134 12424 0 
teezo. 1.575 02055 0.20 ono e176 row 0 
OezZ5 1.969 O2065 OS i269 S2225 29 Ae 0 
foo. 2.362 O23 0:73,0 25362 2) SPATS pNP yi: 0 
toe) 2.756 0.206 O35 DOM RSLS 2555 Pa fis) 3 0 
fea, 3.150 0.3178 0.40 Sa b50 C250 3. Lie 0 
m45 3.543 0.494 0.45 S545 Cay = 30) 32425 0 
we 50° 3.937 0.778 0350 B2o37 0.929 S27 20 0 
feos 4.331 eZ os 0-55 aos ie 3.0 1 4.001 0 
@-60 4.724 2.035 0250 an 724 i. O72 an Zio il 
65 5.118 32433 0.65 5. Vs A sh 4.522 ab 
Ls THO)” Mages ys Be, 6. 000 0.70 Do 2 a) Sov 4.767 iz 
754 15.906 B10.885 OS 52206 7 302d 4.998 2 
Broce  6:299 21.134 0.80 baz o54 Fee. 206 oy ere AWS, 5 
fees 6.693 45.254 O65 62693 m75- 264 aye OU, 4 
0.90 7.087 99.999 C290 7 0C 7a mol. 641 By als) 25) s 
m5 7.480 99.999 0.35 7.480 99.999 Seer al, 7 
moo 7.874 99.999 1.06 oo) ee oe 62005 se) 
mern92. . . SADAIT.: end of Tisting 


5 


S005 
= Vem 
EO 
.144 
a ges 
-245 
poe Al, 
.416 
~o42 
.706 
928 
200 
.564 
039 
655 
-464 
gol 
Poe © 
Oo ae 
.000 


U.S. Army CECOM@MRVECD FLIRIZ Thu, Maye 2s: 32-19 e 200 


output file: SADA1.1 Shore listing 

data file: SADAITI 

command line arguments: -d SADAII -o SADA1 -p MDT 

begin data file listing... 

gen2: sample data file for 2nd generation FLIR with SADA II FPA 


MDTD AT 300 K BACKGROUND TEMPERATURE 


ea MDTD 
O22 0s ese 4o7 05) 27.084 
0.40 19.685 62oG7 
O760 Siege 23 B28 
0.80 9.843 Tole 
O00 12.874 AL AOS 
ie? O 6a >0'2 O87 2 
iy, 20 SB. 624 07672 
1.60 4.921 0.540 
1.80 A374 ol aso 
2.00 be 93 7 0.384 
220 a Ses: 02335 
Ze ope oa Oe296 
2260 B028 0.266 
2260 2612 Ge241 
2206 2625 G22 1. 
3226 2.461 G2 04 
3.40 ve SE) ah. Ae 
260 Paulas) Oealag i s' 
3.680 2072 0.167 
4.00 I Sere, OSs 
4.20 Iss) Or 20 
4.40 1.720 0.142 
4.60 im 712 Oe 36 
4.80 1.640 Or 30 
5.00 i 572 (Oa as 
FLERO2Z ere SODAL ae end of Jistumao 
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APPENDIX B. FIGURES OF COMPARISON OF FLIR92 MODEL SENSOR 
OUTPUTS WITH THE OTHER SENSORS IN WinEOTDA MODEL 
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IFOVx 
(horizontal 


1 oo 0.175 _fo.o005__—fo.o005_ 06 = 5 99.9988 
2_|101_ 0.175 _|o.0005__—o.0005_—p.06 =~ ts 99.9988 
3 l102_ 0.539 [0.0004 fo.0004_—fo.25 0.31 
5 |104 0.4515 (0.0005 0.0005 __—0.25 
6 |105 0.3815 0.0008 _—fo.o008 fo. ft sf 981 2 
7 106 0.2695 0.00066 __fo.00066__ 0.19 0.12 14 99.9 
8 l107_ [0.35 _‘fo.00105 0.00105 .02_ 0.16 ~f.77008 100 
9 {108 [0.1575 0.0005 fo.000374_ 0.05 ot 
13 [112 |0.2485 [0.000307 0.000306_0.225 0.0242. 
16/115 
18 [117 (0.329 0.00095 0.00113 0.04 ~~ os. —S— 
21 120 0.1995 [0.0012 _—fo.0015__—io.067__—io.oase7_—fo6 7 
24 [123 (0.1015 [0.000128 0.000104 0.1  ~—s oon ~—S2 
27 126 (0.175 o.0005_—o.0005_—f.06 =~ 99.9988 


Table B.1 - WinEOTDA and FLIR92 Sensors Physical and 
Performance Parameters 
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Figure B.1 - WinEOTDA and FLIR92 Sensors NETD Comparison 


IFOVx vs Sensor number 


lIFOVx (mrad) 
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Figure B.2 - WinEOTDA and FLIR92 Sensors Horizontal IFOV 
Comparison 
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[FOVy vs Sensor number 
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Figure B.3 - WinEOTDA and FLIR92 Sensors Vertical IFOV 
Comparison 





MRT min vs fs 


MRT min 
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Figure B.4 - WinEOTDA and FLIR92 Sensors Minimum MRT vs. 
Spatial Frequency 
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MRT max vs fs max 
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Figure B.5 - WinEOTDA and FLIR92 Sensors Maximum MRT vs. 
Spatial Frequency 
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APPENDIX C. SEARAD INPUTS 


acne hl llCLlUe 
| [MODEL —“‘“C;CSCOC*™*™C™C*d Mld-Tatittuule summer 
i ietwaan te cicitues 
ATMOSPHERIC PATH between two altitudes 
|_| TEMSCT DETERMINES THE MODE OF EXECUTION |Transmittance mode 
dee = Ee Eee 
WITHOUT MULTIPLE SCATTERING scattering 
TBOUND (K) IS THE BOUNDARY TEMPERATURE 294 
FOR SLANT PATH THAT INTERSECTS THE EARTH 
OR GREYBODY 
2; INPUT FOR CARD 2) OSOC=“CSsS—“sSC—‘“‘“‘<;72S~; DSCC 


aah elon acN, L.VULCN, fESTL, IChDReivycA, 3=Navy maritime 
VIS, WSS, WHH, RAINRT, GNDALT IHAZE extinction, set own 





SeUbGioe Lon BATINGI EON TYPE AND THE Visio 1 Gs 


DEFAULT VISIBILITY RANGE 
ICSTL IS THE AIR MASS CHARACTER (1 TO 10) | 3=open ocean 
. | USED ONLY WITH NAVY MARITIME MODEL naw. Lae 
[ssw sass to seuss Pe RE 
RAIN RATES TO BE USED 


|_| VIS SPECIFIES THE METEOROLOGICAL RANGE | 24.14km=15mi 

eee a 
(AVAILABLE ONLY WHEN IHAZE=3/10) 

ee 


|_| Hi - SPECIFIES THE INITIAL ALTITUDE 
| | H2 SPECIFIES THE FINAL ALTITUDE co 


| | RANGE SPECIFIES THE PATH LENGTH (KM) 0, 
AT THE PARTICULAR GEOGRAPHICAL LOCATION —————, 
|_|xoprrrearrom nu Se USED 
MODIFICATIOIN WILL BE USED 
ee 


INITIAL FREQUENCY (WAVENUMBER CM-1) 1000 
FINAL FREQUENCY (WAVENUMBER CM-1 ) 1255 


‘= = 
¥ DV = = FREQUENCY INCREMENT (OR STEP SIZE) 5 
-. IFWHM = INCREMENTAL FREQUENCY WIDTH AT 10 


Table C.1 - SeaRad Input Parameters 





01. 
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APPENDIX D. SEARAD OUTPUTS 
wkkekek* SEARAD kKkKKK* 


DATE: 05/16/2008 Wane 01:41:29,225 
TRANSMITTANCE MODE 
SINGLE SCATTERING USED 
MARINE AEROSOL MODEL USED 
WIND SPEED = 2.5/7 M/SEC 
WIND SPEED 2.57 M/SEC, 24 HR AVERAGE 
RELATIVE HUMIDITY = 76.11 BEReGENT 
AIRMASS CHARACTER 328 
Vises lhi Ty Za Veer 


SLANT PATH, (Hab TO Hz 


Bi = .010 KM 
H2 = .500 KM 
ANGLE = .000 DEG 
RANGE = .500 KM 
BETA = .000 DEG 
LEN = 0 


FREQUENCY RANGE 


nV = 630 °CM-1 (123 0ee Per OMe TERS) 
rVyZ = 1 25 0.CM—-1. (o2 00s meROME TERS } 
TDW = 5 Cie 

IFWHM = 10 CM-1 

TPLuTER = 0 


SUMMARY OF THE GEOMETRY CALCULATION 


iil = oso KM 

HZ = .500 KM 
ANGLE = Ji 479 DEG 
RANGE = .500 KM 
BETA = .001 DEG 
Pra = 168.521 DEG 
HMIN = 010 KM 
BENDING = .000 DEG 
LEN = 0 


INTEGRATED ABSORPTION 64.47 CM-1 FROM 830010 > 1250 GM 
AVERAGE TRANSMITTANCE = .8465 


#03 







(ian enema | 









10 O21386 


Table D.1 - SeaRad Outputs for Different Ranges 









Transmittance vs Range 
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Figure D.1 - Searad Outputs For Mid-Latitude Summer Scenario 
Ceneitions 
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APPENDIX E. WINEOTDA MODEL DATA ENTRY FORMS 








| wi, WinEGTDA - Targets 


Target [Gunboat v , 


Heading 





= Cancel | 


Figure E.1 - WinEOTDA Target Entry Form 





| ia. Wreainiy. Backgrounds | BEE x] 


Background #t1 | Background #2 | _Backoround #3 _| 
Background Type : [water | 

Depth (fy: [20 | 

Clarity : [Cer + 

nAft 8 =——«‘— 

Slepe: ‘Value sue Direction | 


OK | Cancel | 


Figure E.2 -— WinEOTDA Background Entry Form 





O'S 


ia, WinEOTDA - Time Over Target 


Planning Intervak © 75 © 30 © 66 


TOT - Date: [01/07/1993 Hr: [10 Min:|55 


January 1993 [January ~| [1933 | 





Figure E.3 -—- WinEOTDA Time Over Target (TOT) Entry Form 


Sensor ID:]127 + 


Sensor Height (ft AGL) 
[2000 = 


Scene Complexity 
f* None 
Low 
™ Medium 
™ High 





Cancel | 


Figure E.4 —- WinEOTDA Sensor Entry Form 
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ESR Ses sss es a 
[Time GMT (His} | 10] 47] 12) 13] 14] 15) 16] 17] 18) 13} 20) 21} 22} 23] oo} at 3 


=, (deg F) 51 49 47 #45 45 -45 47 #49 51 55 57 61 64 6 69 71 71 71 
| Dew PT Temp (deg F) 47 45 44 43 43 #43 44 #45 #47 #49 «50 51 53 55 5S? 5? 5S? 57 
Relative Humidity (%) 84 8/f 90 92 93 92 90 SF 84 SO 76 72 69 6 65 G4 64 64 


Agrosol index ===} " 52) 526 2] S22) 2] 2) 82] 22] 2] 2) 2), 92) 2) ai af 2a} 2 








Precipitation Index 
at-abat=\- 


| Wind Speed {kts} 2 EE EE . 15/1 
0} 190 190] 190 ae 


Nested ——eorer aa cra SEER OCeESeS 

Visibility (ri) | 6.2) 2) 62) 6.2) ees ao 6. 
0 
0 





| Wind Direction (deg}__| 130] 130] 130] 130] 130 130 298] 5 190] 90| 190] 13 
Sanday eH 150 150 150 150 150 150 150 150 150 150 150 150 150 150 1 150 150 150 150 15 
ener an Ot) ate Gl td So) Be i SP Po ee i Pe 
LowCloudHt{hft) | Of} Oo] Of Of] 40) 40] 40] 40) 40} 40] 40] 40) 40) 40] 40] 40] 40] 4 
CSSD 
“ Cloud Type 3 

Mid Cloud Ht (hft) 13 va 8 c r F r r a i i rc ta 12 
Mid Cloud Amount (8ths) 1 1 


High Cloud Type one oon eee aS me 
High Cloud Ht (hit) 250) 250] 250] 250] 250/ 250) 250] 250] 250] 250) 250] 250] 250 mates 250] 250 25 


High CloudAmount (Sths}!_2| 2] 2] 2} at 2) 2] ttt} 222 
Time GMT (Hrs) | 10} 11} 12] 13] 14] 15] 16) 17] 18] 13} 20) 21{ 22] 23) oO} 1] 2] 3] 4 






22 SS ee Sessa 
i ee 
Ste Min Téa: [45 Max Tew: a 


| Cancel » | 





Figure E.5 - WinEOTDA Meteorology Data Entry Form 
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APPENDIX F. WINEOTDA MODEL OUTPUTS 


IR EOTDA EXECUTION SUMMARY 


SYSTEM INPUT FILE NAME: C:\PROGRAM FILES\NRL- 
MONTERE Y\WINEOTDA\Data\local\EOTDA\State\S ystemState.dat 


TOT : 07 January 1993 1055 GMT (Z) 
Absolute Humidity : 8.9 (g¢/m**3) 4 km Transmissivity : 0.60 
Sky penis : 243.4 (deg K) an Visibility : 066.3 (kit) 
teenie : 24 deg 15 minN renee : 059 deg 45 min W 
Sensor ID 2 Sensor Ht : 2,000 feet 

View Direction : 0 eee : None 

TARGET INFORMATION 

Target Elevation: 0 feet (MSL) 

Target #1 Target #2 

Target ID : Gunboat Target ID : Gunboat 

Target Heading : 0 Target Heading :0 

Operating State : 1 Operating State :1 

Target Speed :0 Target Speed :0 

BACKGROUND INFORMATION 


General Background Albedo: Ocean 


Background #1 ID: Water 
Background #2 ID: Soil 
Background #3 ID: Vegetation 


mm cm mcm cr cm i re mr em 
mm ni me me a er re es re ee ee ee ee ee ee 





IR EOTDA OUTPUT 

RANGES Target #1 Gunboat Background #1 Water 
View Dir | MRT Detection Range (km) | MDT Detection Range (km) | Lock-on Range 

(deg) | NFOV WFOV | NFOV WFOV | (km) 
woee-en--- |------------|------------- |------------|------------- |--------------- 

000 48.1 18.5 31.8 53.7 0.0 

045 oy es 18.5 303 38.6 0.0 

090 S165 Re) 315 39.5 0.0 

135 51.0 18.5 34.6 38.4 0.0 

180 45.1 18.5 28.8 31.8 0.0 

225 49.1 18.5 30.1 36.9 0.0 

270 49.8 18.5 3h 38.2 0.0 

215 50.4 18.5 31.3 38.0 0.0 
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THERMAL CONTRAST (Delta-T) 


Target #1 Gunboat 


Background #1 Water 





View Dir | MRT Delta T (K) | MDT Delta T (K) | Lock-on Delta-T 
(deg) | NFOV WFOV | NFOV WFOV | 
---------- |-----------|------------- |-----------|------------- |--------------- 
000 Bilis: 20s2 333 15.4 0.0 
045 28.8 25.8 31.0 222 0.0 
090 Ig | 24.2 S12 24.8 0.0 
135 25.9 23.0 28.6 21a 0.0 
180 ZA 18.8 22, 13 0.0 
225 21.0 18.3 27.0 18.0 0.0 
270 21.3 18.5 29.8 20.7 0.0 
eles, 24.4 21.6 B10 20.9 0.0 
TEMPERATURES (K) Target #1 Gunboat Background #1 Water 
View Dir|Bkgd Temp |MRT Temperature (K) !MDT Temperature (K) | Lock-on Temp 
(deg) | (K) INFOV WFOV |INFOV WFOV | (K) 
~--------- |------ |---------|---------- |--------- |---------- |--------------- 
000 Dag ie 304.9 304.6 307.7 289.7 300.0 
045 255 (Va) 3022. 302.2 305.4 296.1 300.0 
090 Pi as 300.5 300.5 305.6 298.7 300.0 
135 2155) 299.4 299.4 3029 e057 300.0 
180 27355 295.3 295.4 297.5 285.9 300.0 
225 DIG) S 294.8 294.9 3018. 22923 300.0 
270 BB 295.0 295.1 304.5 294.9 300.0 
315 2135 298.0 298.0 305.5 294.9 300.0 
RANGES Target #1 Gunboat Background #2 Soil 
View Dir | MRT Detection Range (km) | MDT Detection Range (km) | Lock-on Range 
(deg) | NFOV WFOV | NFOV WFOV | (km) 
---------- |-------------|------------- |-------------|------------- |--------------- 
000 41.0 18.5 26.8 30.1 0.0 
045 42.9 18.5 25:6 32.6 0.0 
090 41.8 18.5 26.4 sles 0.0 
135 40.3 1S 24.7 30.3 0.0 
180 26 16.1 19.8 720 0.0 
225 28.3 17.8 2A, D9) | 0.0 
270 30.1 ikises 26.0 2355 0.0 
315 38.6 18.5 26.4 28.8 0.0 
THERMAL CONTRAST (Delta-T) Target #1 Gunboat Background #2 Soil 
View Dir | MRT Delta T (K) | MDT Delta T (K) | Lock-on Delta-T 
(deg) | NFOV WFOV | NFOV WFOV | 
---------- |-------------|------------- |-------------|------------- |--------------- 
000 10.7 10.5 14.8 6.8 0.0 
045 8.1 8.1 Wes 8.4 0.0 
090 6.4 6.4 13.8 Ve 0.0 
135 a3 SAS, Js’ 6.4 0.0 
180 ‘ 13 6.0 210 0.0 
225 0.8 0.8 DE oe 0.0 
270 0.9 1.0 11.4 2.8 0.0 
315 3.9 3.9 11.4 a3 0.0 


Ie 


mem mm ae me rs em a a a a ee Se ee — ee ee 
SS SS SS SS SS SS SS SS SS SS SS SSS ee ee ce ee ee ee 











TEMPERATURES (K) Target #1 Gunboat Background #2 Soil 
View Dir|Bkgd Temp !MRT Temperature (K) | MDT Temperature (K) | Lock-on Temp 
(deg) : (K) | NFOV WFOV | NFOV WFOV | (K) 
----------|----------- |----------|---------- |----------|---------- |--------------- 
000 294.1 304.8 304.6 308.9 300.9 300.0 
045 294.1 2022 3072 S070") 502.5 300.0 
090 294.1 300.5 300.5 BO7o 301.5 300.0 
[33 294.1 299.4 299.4 305.4 300.5 300.0 
180 294.1 295.4 295.4 300.1 296.1 300.0 
225 294.1 294.9 294.9 B05-2— 296.3 300.0 
270 294.1 295.0 295.1 305.5 296.9 300.0 
By) 294.1 298.0 oe 305.5 299.4 300.0 
RANGES Target #1 Gunboat Background #3 ee con 
View Dir | MRT Detection Range (km) | MDT Detection Range (km) | Lock-on Range 
(deg) | NFOV WFOV | NFOV WFOV | (km) 
--+------- |-------------|------------- |-------------|------------- |--------------- 
000 34.1 18.5 25-0 Zo-8 0.0 
045 31.3 18.5 JE) 21.9 0.0 
090 24.7 15.9 23.4 255 0.0 
135 320 18.5 22.5 25.6 0.0 
180 36.5 18.5 24.9 e735 0.0 
22) 41.0 18.5 29.4 BI.) 0.0 
270 41.2 18.5 PIE 7: 31.1 0.0 
ole 36.9 18.5 26.8 Zio 0.0 
THERMAL CONTRAST (Delta-T) Target #1 rearnnerrs Background #3 eecoean 
View Dir | MRT Delta T (K) | MDT Delta T (K) | Lock-on Delta-T 
(deg) | NFOV WFOV | NFOV WFOV | 
---------- |-------------|------------- |-------------|-------------_|--------------- 
000 3.9 ot 9.1 35 0.0 
045 les 13 T6 2.0 0.0 
090 -0.4 -0.4 9.0 -3.] 0.0 
135 -1.5 -1.5 -7.4 -3.3 0.0 
180 -5.5 -5.5 -7.6 -4.1 0.0 
O25 -6.0 -6.0 -7.6 -6.6 0.0 
270 -5.9 -5.8 -8.4 -6.6 0.0 
BS -2.9 -2.9 -7.6 -4.1 0.0 
TEMPERATURES (K) Target #1 Gunboat Background #3 yeccuicn 
View Dir|Bkgd Temp |MRT Temperature (K) | MDT Temperature (K) | Lock-on Temp 
(deg) | (K) | NFOV WFOV | NFOV WFOV | (K) 
~--------- |----------- |----------|---------- |----------|---------- |--------------- 
000 300.9 304.8 304.6 310.0 304.2 300.0 
045 300.9 302.20 33022 303.5 aeo0Z.9 300.0 
090 300.9 300.5 300.5 309.9 297.8 300.0 
135 300.9 299.4 299.4 £935 297.5 300.0 
180 300.9 295.3 295.4 293.3 296.8 300.0 
E25 300.9 294.8 294.9 293.3 294.2 300.0 
270 300.9 295.0 295.1 292.5 294.3 300.0 
515 300.9 298.0 298.0 293.3 296.8 300.0 


jieae 






a. WinEQTDA OUTPUT | 
File Target Backgrounds Detection Range 








— Gun_Wate_N 






IR EOTDA EXECUTION OUTPUT 





Gunboat Gunboat 
Range {kft) Range (kft) 






View Direction 


106 1979 1106 
1266 1690 1266 
1236 1690 1296 


126.0 one 126.0 __ 







1044 148.1 9104.4 
121.0 161.0 121.0 
125.3 1635 1253 


ot pepeet onan Y hee 


Rieig. 1653) 14n 









0.9 -=-> No Value Computed 
-1.90 --> No Solution Possible 





Figure F.2 - WinEOTDA Tabular Output 


2 


PR wink OTOA [eT] 
Fie Edt input Bun Quip Help 


Olea) S| Hla} slesl) SO] 2] =| Of B| 





a ae P+ ee ea pty egy 





Met Summary ~ 

Temp: 76F RH: 40% ~ 
Wind 120 5.49 kts 

Vis: 14.3 (wi) 
Scattered 






Sensor 





IR 8127 
Time Over Target 2000 ft 0 deg 
Exe Tene: 
01/07/1993 1055 GMT(Z) 
Backgrounds 
ft Ale Targets Gunhoat 
Trea 6 (ee 
; ¥ Ge fect 


x wide: ee 2 f veg aie * Sawer 
es Re en 


-s 


Target Location: 24 15 N 5945 W 





. = : se 
ANAM ae et RA Ml i OREN LN RA RR PN ABET RNAS PARA NRA TAMARA ARMERRR Gtr Te te line Re Rn pS awe 


Status: IR EOTDA Ran Successfully, 


tatew os tote Aenean ame Smee! wee came re at ete Se ene eee een etree tre mnen eeeee | fee eneen tM tmn met Tt reat Mt le amet \iteeervwere vere tentrer sett] Lhe ee et am wees pal tts ale 


Figure F.3 - WinEOTDA Main Screen Output 
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APPENDIX G. ACQUIRE MODEL INPUTS 


Pequice Gaza file for freent view Of gunboat ZO000ftE at 25km 


>sensor_lookup 


data_file_name Sada.dat 

sensor_id gen2 

performance_mode MRT 
>target 

characteristic_size 9.78 

target_signature Tico 


>cycle_criteria 


detection_n50 0.75 
detection_n50 0.75 
ellaissirticativon n50 es: 
recognition_n50 o20 
identification_n50 6.0 


>band-averaged_atmosphere 
#_ points: 20 Ke aes 6s eee 

.000e+00 
.000e-01 
.500e-01 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+01 
.200e+01 
.400e+01 
.600e+01 
.800e+01 
.000e+01 
.200e+01 
.400e+01 
.600e+01 
.800e+01 


NONNNNPRPPRPRPRPOoOO U1 B® WN Fe ul © 


>end 


le} 


MRT_MDT_MRC_OR_MDC 


meters 
degrees _C 


wfov 
nfov 
igvicionts 
Neo 
iQGlons 


transmittance 


SIDRPRPPNWMHAIWOPPRPN WA PHO ~1)~) OF 


.000e+00 
.465e-01 
.960e-01 
.510e-01 
.049e-01 
.942e-01 
.070e-01 
.371e-01 
.804e-01 
.959e-01 
.386e-01 
. 860e-02 
.090e-02 
.110e-02 
./10e-02 
. 700e-02 
.970e-02 
.440e-02 
.060e-02 
.800e-03 
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APPENDIX H. CRITICAL DIMENSION ANALYSIS OF ACQUIRE TARGET 
MODEL 


Ceeom [SRE | Gee bce negnian fae 
3 500 ft 500 ft (Front View Front View 500ft View Side View 500ft 


Table G.1. -— Critical Dimension Analysis for 500ft Sensor 
Height 





























ah] 


rane [Sat | mo alsa el Ss 
: 2000 ft 2000 ft |Front View Front View 2000ft| Side View| Side View 2000ft 
3000 
5000 
104.6794 


30000 609.6 0.020321 (94.00426 |9.69558 373.3466 |19.32218 
Table G.2 - Critical Dimension Analysis for 2000ft Sensor 
Height 


#18 






: 4000 ft | 4000ft (Front View Front View 4000ft|Side View| Side View 4000ft 
COP se) 
24000 [1219.2 _fo.050822 [106.2442 |10.30748 __—__—f885.2834 


Table G.3 - Critical Dimension Analysis for 4000£ft Sensor 
Height 

























daks® 


Critical Dimension (CD) (m) 





0 4000 8000 12000 16000 20000 24000 28000 32000 > 
Range (m) 


—<CD Side View 2000ft —*-CD Front View 4000ft —e— CD Side View 4000it | 


rt a A ER Sl a NE | 


4000ft Sensor Heights 
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APPENDIX I. ACQUIRE MODEL OUTPUTS 


mun #1 
U.S. Army CECOM RDEC NVESD 
ACQUIRE version 1 


sadallfa 
-d sadalIlfa 


data file: 
command line: 


(Mavers 0 1995) 


TARGET DISCRIMINATION RANGE PERFORMANCE 


begin parameter listing... 


Acquire data file for front view of gunboat 500ft at 25km 


>sensor_lookup 
data_file_name 
sensor_id 
performance_mode 


>target 
characteristic_size 
target_signature 


>cycle_criteria 
detection_n50 
detection_n50 
classification_n50 
TeCOOnM. -10l. 1oU 
fdentification n50 


>band-averaged_atmosphere 
# points: 20 


>end 


emer parameter listing. = 


NONNNNRPRPRFRRF PRP ODUM HhWNFH TWO 


Sada.dat 


gen2 
MRT 


939 
tage 


eo e©« ee e© © © @ @ 


.000e+00 
.000e-01 
-500e-01 
-000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+01 
-200e+01 
.400e+01 
.600e+01 
.800e+01 
.000e+01 
.200e+01 
.400e+01 
.600e+01 
.800e+01 


YPRPRPNWUDWOFPENWA BRAY) DOF 


i 2h 


MRT_MDT_MRC_OR_MDC 


neocons 


degrees_C 


wfov 
NEOV 
TWLOV 
nNLov 
nEov 


transmittance 
.000e+00 
-465e-01 
.960e-01 
.510e-01 
.049e-01 
.942e-01 
-070e-01 
.371e-01 
.804e-01 
-959e-01 
.386e-01 
-860e-02 
.090e-02 
.110e-02 
.710e-02 
.700e-02 
.970e-02 
-440e-02 
.060e-02 
.800e-03 


Wed May 31 21:33:59 2000 


MESSAGES 


2D MRTD from lookup table 


cyles/mrad MRTD 

8.300e-001 6.500e-002 
1.120e+000 8.500e-002 
1.424e+000 1.100e-001 
1.767e+000 1.440e-001 
2.112e+000 1.880e-001 
2.458e+000 2.450e-001 
2.794e+000 3.190e-001 
3.116e+000 4.160e-001 
3.425e+000 5.420e-001 
3.720e+000 7.060e-001 
4.001e+000 9.210e-001 
4.268e+000 1.200e+000 
4.522e+000 1.564e+000 
4.767e+000 2.039e+000 
4.998e+000 2.658e+000 
5.219e+000 3.464e+000 
5.429e+000 4.515e+000 
5.629e+000 5.886e+000 
5 .820e+000 7.672e+000 
6.003e+000_ 1.000e+001 


Sky-to-ground ratio defaulted to 1 for thermal systems. 
Last range for input atmospheric transmittance data is less 
than maximum range. 
Loadline does not intersect MRTD/MRC above range 27.4 km, extend 
curve to lower frequencies and temperatures. 
SENSOR 
gen2 from Sada.dat 
TARGET 
characteristic dimension: 9.39 meters 
inherent signature: 7.71 degrees C 


OBSERVER ENSEMBLE PERFORMANCE 
RANGE GIVEN PROBABILITY... 


WFOV NFOV NFOV NFOV NFOV 
prob N50=0.75 N5SO0=0275 Nsg=1.50 N50=3 .00 N50=6.00 
O35 6267 km ie Zoe en oe 76 kn Je O wen 3.94 km 
OB Ede Bae 18.36 Zee Scr 4.51 
0285 10266 15:6 LEC ae S267 4.94 
O2s0 ile 7 os 14.42 9a23 5230 
Ona 5 dato 2022S ta. 99 ye Gil 5.64 
O70 PR 208 7 ops 2 1O0ea 52.96 
O65 eae 2a ise 1EO (3) 1 6.27 
0-60 13524 21.49 16.45 ites 6255 
ers, eT 2ieee 6.91 Ia 6289 
O50 14.18 2220 ee 2 dealt eae 
0.45 14.64 22003 17.62 1 3 SS, aad, 
0.40 15.13 PPR LE, 1629 12266 (5 
0235 iS o5 PX TENG, ier 6 1326 6235 
0.30 I BE: 2365 19'°.33 13m 9S 8.83 
O25 Tie se2 24.30 Ne he. 14.61 9439 
0220 17.64 24.81 20e 50 1S-3e 10.05 
ORs: 18.54 25o47 2 dees BS cil IO = 
O10 1S. 72 263) 22.46 7260 12-05 
0.05 Ze so 0.00 24.07 desl 14.29 


end of run 1 from sadalIIfa 
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Figure H.1. Range Performance Gunboat Front View at 500 ft 
Sensor 


25 


pan bh ole ab Wed May Si921:23:47 2000 
U.S. Army CECOM@BREESG NVESD 
ACQUIRE version 1 (May 30 1995) 


data file: sadalIIf 

command line: -d sadalIIf 

TARGET DISCRIMINATION RANGE PERFORMANCE 

begin parameter listing... 

Acquire data file for front view of gunboat 2000ft at 25km 


>sensor_lookup 


data_file_name Sada.dat --- 

sensor_id gen2 --- 

performance_mode MRT MRT_MDT_MRC_OR_MDC 
>target 

characteristic_size 9.78 meters 

target_signature feel degrees_C 


>cycle_criteria 
detection_n50 G27 wfov 
detection_n50 8 gear nfov 
Glassificatien.n50 dss niov 
wecognitionens 0? 3.0 nfov 
poentification n50 60 nfov 


>band-averaged_atmosphere 


#_points: 20 | ee eee transmittance 
0.000e+00 1.000e+00 
5.000e-01 8.465e-01 
7.500e-01 7.960e-01 
1.000e+00 7.510e-01 
2.000e+00 6.049e-01 
3.000e+00 4.942e-01 
4.000e+00 4.070e-01 
5.000e+00 3.371le-01 
6 .000e+00 2.804e-01 
8.000e+00 1.959e-01 
1.000e+01 1.386e-01 
1.200e+01 9.860e-02 
1.400e+01 7.090e-02 
1.600e+01 5.110e-02 
1.800e+01 3.710e-02 
2.000e+01 2.700e-02 
2.200e+01 1.970e-02 
2.400e+01 1.440e-02 
2.600e+01 1.060e-02 
2.800e+01 7.800e-03 


>end 


end parameter listing... 


124 


MESSAGES 


ZO=MRID Exom Ieekup stable 


cyles/mrad 
.300e-001 
.120e+000 
.424e+000 
. 767e+000 
.112e+000 
-458e+000 
.794e+000 
.116e+000 
.425e+000 
.720e+000 
.001e+000 
.268e+000 
.522e+000 
.767e+000 
.998e+000 
.219e+000 
-429e+000 
-629e+000 
. 820e+000 
.003e+000 


NNN NN FP PPP SP WWW Dd DN NP EF PF OO 


PAU PWNNYPPRP WO DUO BP WNRPP BOD 


MRTD 


.500e-002 
.500e-002 
.100e-001 
-440e-001 
.880e-001 
-450e-001 
.190e-001 
.160e-001 
-420e-001 
.060e-001 
.210e-001 
.200e+000 
.564e+000 
.039e+000 
.658e+000 
.464e+000 
.515e+000 
.886e+000 
.672e+000 
.000e+001 


Sky-to-ground ratio defaulted to 1 for thermal systems. 
Last range for input atmospheric transmittance data is less 
than maximum range. 
Loadline does not intersect MRTD/MRC above range 27.4 km, 
curve to lower frequencies and temperatures. 


SENSOR 


gen2 from Sada.dat 


TARGET 


characteristic dimension: 
inherent signature: 7.71 degrees C 


OBSERVER ENSEMBLE PERFORMANCE 
RANGE GIVEN PROBABILITY... 


WFOV 
BEOD NSO=0.75 
0295 92145 len 
0.90 10220 
Sac Oo 7 
0.80 57, 
oS 5 Eel 
6.70 2365 
0.65 eS ceO 
0.60 S50), 
Oe 55 14.04 
6250 14.49 
o245 14.96 
0.40 15 47 
eo. 35 cS. 99 
9.30 £6.55 
G25 dey eee 
6.20 17.94 
6.15 e286 
o- 10 ZO 0 
o-05 Ziao> 


NFOV 
N50=0.75 


end of run 1 from sadalII£ 
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Figure H.2. Range Performance Gunboat Front View at 2000 ft 
sensor 


LATE 


Cun 
Wes. Army CECOM RDEC NVESD 
ACQUIRE version 1 (May 30 1995) 


data file: sadaIIfb 
command line: -d sadaIIfb 


TARGET DISCRIMINATION RANGE PERFORMANCE 


begin parameter listing... 


Acquire data file for front view of gunboat 4000ft at 25km 


>sensor_lookup 


data_file name Sada.dat 

sensor_id gen2 

per formance_mode MRT 
>target 

characteristic_size Ore? 

target_signature Fert 


>cycle_criteria 
detection_n50 07> 
detection_n50 CaS 
classification_n50 1 
recogiit LOMens]® ie 
Laentificaticon= m5” 6 


>band-averaged_atmosphere 


# points: 20 kn. . ee 


.000e+00 
.000e-01 
.500e-01 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+00 
.000e+01 
.200e+01 
-400e+01 
.600e+01 
.800e+01 
-000e+01 
.200e+01 
.400e+01 
.600e+01 
.800e+01 


NONMNMNFFRP RPP ODUM RP WNP DOO 


>end 


end parameter listing... 


...transmittance 
.000e+00 
.465e-01 
-960e-01 
.510e-01 
.049e-01 
-942e-01 
.070e-01 
.371e-01 
.804e-01 
-959e-01 
.386e-01 
.860e-02 
.090e-02 
-110e-02 
.710e-02 
-700e-02 
.970e-02 
-440e-02 
.060e-02 
.800e-03 


JP RPRPNWUDWOPPRPNWA RB OD YD~I OP 


Lag 


MRT_MDT_MRC_OR_MDC 


meters 


degrees_C 


wfov 
nfov 
MEOV 
nfov 
mWLOV 


Wed May 31 21:34:14 2000 


MESSAGES 


2D MRT Deerom Lookup table 


cyles/mrad MRTD 

8.300e-001 6.500e-002 
1.120e+000 8.500e-002 
1.424e+000 1.100e-001 
1.767e+000 1.440e-001 
2.112e+000 1.880e-001 
2.458e+000 2.450e-001 
2.794e+000 3.190e-001 
3.116e+000 4.160e-001 
3.425e+000 5.420e-001 
3.720e+000 7.060e-001 
4.001e+000 9.210e-001 
4.268e+000 1.200e+000 
4.522e+000 1.564e+000 
4.767e+000 2.039e+000 
4.998e+000 2.658e+000 
5.219e+000 3.464e+000 
5.429e+000 4.515e+000 
5.629e+000 5.886e+000 
5.820e+000 7.672e+000 
6.003e+000 1.000e+001 


Sky-to-ground ratio defaulted to 1 for thermal systems. 
Last range for input atmospheric transmittance data is less 
than maximum range. 
Loadline does not intersect MRTD/MRC above range 27.4 km, 
curve to lower frequencies and temperatures. 

SENSOR 
gen2 from Sada.dat 

TARGET 
characteristic dimension: 10.27 meters 
inherent signature: 7.71 degrees C 


extend 


OBSERVER ENSEMBLE PERFORMANCE 
RANGE GIVEN PROBABILITY... 


WFOV NFOV NFOV NFOV NFOV 
prob N50=0.75 N50=0.75 N50=1.50 N50=3 .00 N50=6.00 
O235 9.48 km 17-89 ion 12.46 km 7 aos ikm 4.26 km 
Gao foe 56 139206 eG Sso7 4.88 
6285 aS 2 Ic be 7/3) 14.48 Jaa See 
0.80 Vi. 94 Z0eoc Ieye Ee, eae: DelZ 
0275 Sleames al 20287 Wo 2 10.34 6.08 
0270 i200 PLA IAS | ioe L 10283 6.42 
O65 S42 Zl 69 E6769 ee 6.74 
0260 1307 220g be 0 eg Lee 7206 
0255 14.42 22.46 io dh 3 US 7.40 
0250 14.87 22.80 18.04 i261 a gS 
0.45 ES.43'6 Zo 5 Te 250 13207 See 
0.40 Te 6 25.53 Lee? 1350 8250 
O23:5 16236 ZS 19.44 Seba 8.94 
0230 16.95 24.33 Ug2 26 14.67 9.44 
O25 L729 24.78 20S IS oa, IO ATON®, 
Op. 8 oe Sal 25250 2s 16.08 IOS 720, 
O25 JS 18, 25.94 2s 97 ec dis 3 
O10 20636 PAS a IES PAP SM ioe e7 12.84 
OOS: 22506 0.00 24.55 20525 14.99 


end Of yun 1 from sadaleco 


As 
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Figure H.3. Range Performance Gunboat Front View at 4000 ft 
Sensor 


ae 


run #1 Wed May 31 21:34:23 2000 
U.S. Army CECOMS@REEC NVESD 


ACQUIRE version 1 (May 30 1995) 


data file: sadalIIsa 

command line: -d sadalIIsa 

TARGET DISCRIMINATION RANGE PERFORMANCE 

begin parameter listing... 

Acquire data file for side view of gunboat S00ft at 25km 


>sensor_lookup 


data_file_name Sada.dat --- 
sensor_id gen2 --- 
per formance_mode MRT MRT_MDT_MRC_OR_MDC 
>target 
characteristic_size iS ygik, 7) meters 
target_signature Youn degrees_C 
>cycle_ criteria 
detection_n50 O75 wftov 
detection_n50 O2TS nfov 
classification. nov aS niov 
recognition_n50 Bi) nfov 
PAentitication n50 Oru nNreov. 
>band-averaged_atmosphere 
fepolnes. 20 HM... .... . aes transmittance 
0.000e+00 1.000e+00 
5.000e-01 8.465e-01 
7.500e-01 7.960e-01 
1.000e+00 7.510e-01 
2.000e+00 6.049e-01 
3.000e+00 4.942e-01 
4.000e+00 4.070e-01 
5.000e+00 3.371le-01 
6.000e+00 2.804e-01 
8.000e+00 1.959e-01 
1.000e+01 1.386e-01 
1.200e+01 9.860e-02 
1.400e+01 7.090e-02 
1.600e+01 5.110e-02 
1.800e+01 3.710e-02 
2.000e+01 2.700e-02 
2.200e+01 1.970e-02 
2.400e+01 1.440e-02 
2.600e+01 1.060e-02 
2.800e+01 7.800e-03 


>end 


end parameter listing... 


130 


MESSAGES 


2D MRTD from lookup table 
cyles/mrad 


AuMNUNUM HNP PP BP BWWWNNNHPE Bo 


.300e-001 
.120e+000 
.424e+000 
.767e+000 
.112e+000 
-458e+000 
.794e+000 
.116e+000 
.425e+000 
.720e+000 
.001e+000 
.268e+000 
.522e+000 
.767e+000 
.998e+000 
.219e+000 
-429e+000 
-629e+000 
.820e+000 
.003e+000 


PHYO PRPWDNDNDP PO DU BWDP Bo 


MRTD 
.500e-002 
.500e-002 
.100e-001 
-440e-001 
.880e-001 
-450e-001 
.190e-001 
.160e-001 
-420e-001 
.060e-001 
.210e-001 
.200e+000 
.564e+000 
.039e+000 
-658e+000 
-464e+000 
-515e+000 
.886e+000 
.672e+000 
.000e+001 


Sky-to-ground ratio defaulted to 1 for thermal systems. 
Last range for input atmospheric transmittance data is less 
than maximum range. 
Loadline does not intersect MRTD/MRC above range 27.4 km, 
curve to lower frequencies and temperatures. 

SENSOR 
gen2 from Sada.dat 

TARGET 
characteristic dimension: 19.17 meters 
inherent signature: 7.71 degrees C 


extend 


OBSERVER ENSEMBLE PERFORMANCE 
RANGE GIVEN PROBABILITY... 


WFOV NFOV NFOV NFOV NFOV 
prob NS50=0.75 N50=0.75 N50=1.50 N50=3.00 NS0=6 . 00 
O95 14.18 km 22.27 ska eo 9 kin 1.9) km aoe. Ie 
290 Tao AGH ike) ree oS 3 < LO 8.14 
0235 16220 232.50 oe 30 13.94 8.80 
0.80 16-69 24.28 9 88 14.58 9237 
eS 17.44 24.66 20.40 L515 9235 
ee 7 0 ia PRs) SOP 20.8) 13669 10.31 
ve 65 18 337 25405 ele 2a LGM LOR a7 
07260 £62835 25.66 2 03 16M61 eyo. 
Oaks Le ee) Paes wees ie Oe L7m0's Tie 62 
ge 50 tec anche: 26n26 22.40 35 Ie 0) 
0.45 20.03 26.56 2273 od 12.54 
0.40 20228 26.88 Zone 18.45 e072 
0) eS Bes PAU Io, PPS PA® Zoo heey 1254 
oe 30 21.44 0.00 23.95 19.47 14.14 
va2D PAL TS, 0.00 24.40 ZU. 05 14.77 
Ue 20 2260 O2.00 24.93 20.74 Als ee fe) 
G15 2351 0.00 25256 PAM ah 16.48 
10 24.26 0.00 26241 PP Mres a | ey ae 
O05 2575 0700 0200 24.18 i975 


end of run 


1 from sadalIsa 


ioe 
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Figure H.4. Range Performance Gunboat Side View at 500 ft 
sensor 


toe 


run FL Wed May 31 21:27:14 2000 
U.S. Army CECOM RDEC NVESD 


ACQUIRE version 1 (May 30 1995) 


data file: sadaIIs 

command line: -d sadaIIs 

TARGET DISCRIMINATION RANGE PERFORMANCE 

begin parameter listing... 

Acquire data file for side view of gunboat 2000ft at 25km 


>sensor_lookup 


data_file_name Sada.dat --- 
sensor_id gen2 --- 
performance_mode MRT MRT_MDT_MRC_OR_MDC 
>target 
characteristic_size 19.36 meters 
target_signature Tere degrees_C 
>cycle_criteria 
detection_n50 O75 wfov 
detection_n50 O75 nfov 
Classification nmo0 es nfov 
recognition_n50 By 0 muhow 
identification_n50 6.0 plage nts 
>band-averaged_atmosphere 
#_ points: 20 Ke. eae ee transmittance 
0.000e+00 1.000e+00 
5.000e-01 8.465e-01 
7.500e-01 7.960e-01 
1.000e+00 7.510e-01 
2.000e+00 6.049e-01 
3.000e+00 4.942e-01 
4.000e+00 4.070e-01 
5.000e+00 3.371le-01 
6.000e+00 2.804e-01 
8.000e+00 1.959e-01 
1.000e+01 1.386e-01 
1.200e+01 9.860e-02 
1.400e+01 7.090e-02 
1.600e+01 5.110e-02 
1.800e+01 3.710e-02 
2.000e+01 2.700e-02 
2.200e+01 1.970e-02 
2.400e+01 1.440e-02 
2.600e+01 1.060e-02 
2.800e+01 7.800e-03 


>end 


end parameter listing... 


L328 


MESSAGES 


2.) MRTDStrom lookup table 


cyles/mrad 
.300e-001 
.120e+000 
.424e+000 
. 767e+000 
.112e+000 
.458e+000 
.794e+000 
.116e+000 
.425e+000 
.720e+000 
.001e+000 
.268e+000 
.522e+000 
.767e+000 
.998e+000 
.219e+000 
.429e+000 
.629e+000 
.820e+000 
.003e+000 


AWN bP PPP BWW WYONYNNPP Bo 


MRTD 
6.500e-002 
8.500e-002 


1.100e-001 


1.440e-001 
1.880e-001 
2.450e-001 
3.190e-001 
4.160e-001 
5.420e-001 
7.060e-001 
9.210e-001 
1.200e+000 
1.564e+000 
2.039e+000 
2.658e+000 
3.464e+000 
4.515e+000 
5.886e+000 
7.672e+000 
1.000e+001 


Sky-to-ground ratio defaulted to 1 for thermal systems. 
Last range for input atmospheric transmittance data is less 
than maximum range. 
Loadline does not intersect MRTD/MRC above range 27.4 km, 
curve to lower frequencies and temperatures. 


SENSOR 


gen2 from Sada.dat 


TARGET 


characteristic dimension: 
inherent signature: 


OBSERVER ENSEMBLE PERFORMANCE 
RANGE GIVEN PROBABILITY... 


WFOV NFOV 
pzeb N50=0.75 N50=0.75 
0.95 14.26 km 22.35 mM 
0290 15.48 23224 
0.85 6.25 23.250 
0.80 16.94 24.53 
0275 yep 24.71 
0570 7 3 Se: Zo 07 
Oves 18.45 25.41 
o) slo 18.91 Zon 73 
0.355 We 2 26.02 
O50 LOZ 26.31 
0.45 ZUraiS PAS (aval 
0.40 20.60 26.92 
0.35 208 vd as 
O30 aA Sal 0-00 
On 25 22.06 0.00 
O20 222.66 0-00 
e.15 23236 0.00 
€.10 24232 0200 
G.05 25.50 0.00 


end of run 


1 from Sadalls 
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19.36 meters 
7.71 degrees C 


NFOV 
N50=3.00 


Side View at 2000ft Altitude 
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Figure H.5. Range Performance Gunboat Side View at 2000 ft 
sensor 


seas 


Tare |: Wed May 31 21:34:32 2000 
U.S. Army CECCM BPEC NVESP 


ACQUIRE version 1 (May 30 1995) 


data file: sadaIIsb 

command line: -d sadaIIsb 

TARGET DISCRIMINATION RANGE PERFORMANCE 

begin parameter listing... 

Acquire data file for side view of gunboat 4000ft at 25km 


>sensor_lookup 


data_file name Sada.dat --- 
sensor_id gen2 --- 
performance_mode MRT MRT_MDT_MRC_OR_MDC 
>target 
characteristic_size Aerts. iL meters 
target_signature Teo7 degrees_C 
>cycle_ criteria 
detection_n50 O25 wfov 
detection_n50 O75 neo 
elassification n50 1S nLeov 
recognition_n50 ae O) niLov 
Peentatication_n50 6.0 rneOov 
>band-averaged_atmosphere 
#_ points: 20 1.91) oa ace transmittance 
0.000e+00 1.000e+00 
5.000e-01 8.465e-01 
7.500e-01 7.960e-01 
1.000e+00 7.510e-01 
2.000e+00 6.049e-01 
3.000e+00 4.942e-01 
4.000e+00 4.070e-01 
5.000e+00 3.371e-01 
6.000e+00 2.804e-01 
8.000e+00 1.959e-01 
1.000e+01 1.386e-01 
1.200e+01 9.860e-02 
1.400e+01 7.090e-02 
1.600e+01 5.110e-02 
1.800e+01 3.710e-02 
2.000e+01 2.700e-02 
2.200e+01 1.970e-02 
2.400e+01 1.440e-02 
2.600e+01 1.060e-02 
2.800e+01 7.800e-03 


>end 


end parameter listing... 


13.6 


MESSAGES 


2D MRTD from lookup table 


cyles/mrad MRTD 

8.300e-001 6.500e-002 
1.120e+000 8.500e-002 
1.424e+000 1.100e-001 
1.767e+000 1.440e-001 
2.112e+000 1.880e-001 
2.458e+000 2.450e-001 
2.794e+000 3.190e-001 
3.116e+000 4.160e-001 
3.425e+000 5.420e-001 
3.720e+000 7.060e-001 
4.001e+000 9.210e-001 
4.268e+000 1.200e+000 
4.522e+000 1.564e+000 
4.767e+000 2.039e+000 
4.998e+000 2.658e+000 
5.219e+000 3.464e+000 
5 .429e+000 4.515e+000 
5.629e+000 5.886e+000 
5 .820e+000 7.672e+000 
6.003e+000 1.000e+001 


Sky-to-ground ratio defaulted to 1 for thermal systems. 
Last range for input atmospheric transmittance data is less’ 
than maximum range. 


Loadline does not intersect MRTD/MRC above range 27.4 km, extend 
curve to lower frequencies and temperatures. 
SENSOR 
gen2 from Sada.dat 
TARGET 
characteristic dimension: 19.61 meters 
inherent signature: 7.71 degrees C 
OBSERVER ENSEMBLE PERFORMANCE 
RANGE GIVEN PROBABILITY... 

WFOV NFOV NFOV NFOV NFOV 
prob NS5O0=0275 NSO=02ys N50=17 50 NS0=300 N50=6.00 
eee e RK He ee ee een ee en eee nee ee ee eer nee eee reer eee He ee UO". 
14.35 km 22.42 km 17.555 san 12.09@ta0 f230° km 
0.90 15258 Zool be. 7 1 L3 226 Sage 
0.85 16.38 25-94 19.46 14.12 8.95 
0.80 0S 24.40 204 L4G 9.52 
Oe 75 1772-60 24.78 2057 LS yoe LO 20m 
or. 0 18.08 Palsy. lel 20.99 Seo Oras 
105 Le.55 25.48 2139 16s PO; os 
0.60 19h OU 25.60 Zo 16.79 Ti oG 
O55 19.40 26.09 22217 17 228 11.80 
0250 19261 26.36 22.53 7 270 £2225 
O45 20.24 26.68 22.88 18.15 Lie al 
0.40 20.69 26.99 2325 18.62 13.19 
O35 2 les? Ve re oil 23-65 iS Pe AL) de ane 0 
Be 30 212538 0.00 24.08 9263 geil 
0225 222155 0.00 24.53 20221 14.95 
D. 20 22. 74 0.00 2505 20.88 LS 
O15 23.44 0.00 2572-7 1. Zoo i665 
0210 24.38 0.00 26555 22./i 7.92 
O05 25-286 0.00 0200 24.30 19.91 


end of run 1 from sadalIsb 


dese) 


Side View at 4000ft Altitude 
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Figure H.6. Range Performance Gunboat Side View at 4000 ft 
sensor 


Leis) 


APPENDIX J. COMPARISON OF BEER’S LAW AND SEARAD OUTPUTS 


BEER ’S LAW APPROXIMATION FOR MIDLATITUDE 
SUMMER OUTPUTS OF SEARAD 


-yR 


The Beer's Law gives transmissivity (:) as: t= e 


ee 
Ip 
The 4km transmissivity calculated by SeaRad is 14 := 0.407 


The atmospheric extinction coefficient (1) for 4km range (R) R:=4 canbe computed as: 


~In(t i) 
R 





te t* In(t4) = -u-R “c= u = 0.225 


Then fory=0.225 , the other transmissivities for different ranges can be found by using the 
Beer's Law approximation. (Ranges, r,, are 0,0.5,0.75,1...6,8...28 km) 





I:= 20 i:=0..I-1 
0 
0.5 0.8465 
0.75 0.796 
0.751 
2 0.6049 
3 0.4942 
4 0.407 
5 0.3371 
6 0.2804 
g 09 0.1959 
a ns PEE ete 
12 0.0986 
14 0.0709 
16 0.0511 
18 0.0371 
20 0.027 
22 0.0197 
24 0.0144 
26 0.0106 
28 Wee | 0.0078. 


IL She, 


Transmissivity 





Comparison of Transmissivities Calculated Using Beer’s 
Law Approximation and SeaRad (Common 4km 
Transmissivity) 


—e— Transmissivities 
calculated using 
Beer’s law 
approximation 

—s— Transmissivities 
calculated using 
SeaRad 






Laks Ue HERELE 
SURE RURSERA Kas 3s <2 
0 10 20 30 

Range (km) 
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APPENDIX K. COMPARISON TABLES OF WINEOTDA AND TAWS OUTPUTS 


A. TAWS OVERVIEW 


The Target Acquisition Weather Software (TAWS) 
predicts the performance of air-to-ground electro-optical 
weapon and navigation systems. The underlying algorithms 
eceumcent ical Eo sthose Cf HOTDS VV. 3.1 (ane WinkOTDA, 
although some programming errors have been corrected. 
Performance is expressed primarily in terms of maximum 
detection or lock-on range. Results are displayed in 
graphic and tabular formats. The program is available 
PeaOugh NRL or throwgm AFRL- 


TAWS supports systems in three regions of the 


Spectrum: Infrared (3-5 micrometers; 8-12 micrometers); 
versible Cea ~ 0.9 micrometers) ; and Laser (Gi OG 
micrometers). The Visible includes both television (TV) and 


Night Vision Goggles (NVG) systems. 
TAWS 1S designed to provide several types of analyses: 
e Illumination Analysis: involves the computation of 
solar and lunar ephemeris information for a 
specified location. A mission planner, for example, 


might be interested in an illumination analysis to 
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determine the time of sunset for a particular 


MLSS MmmaGdahe. andwLocakwens 


e Single Point-Based Analysis: involves detailed 
performance predictions for a particular location. A 
mission planner, for example, might be interested in 
a point-based analysis to predict detection range 
for a particularly important target as a function of 


amen 


e Multiple Map-Based Analysis: involves detailed 
performance predictions for locations along a 
mission route. A mission planner, for example, might 
be interested in a map-based analysis to predict 
detection range for a series of key locations as a 
mimet lon Of Eimes 


TAWS runs on a PC under Microsoft Windows 95/NT/98. 


B. COMPARISON OF TAWS AND WinEOTDA OUTPUTS 


TAWS was run with the same scenario input parameters 
used in WinEOTDA to observe the differences in delta T 
calculations and detection ranges. The delta T outputs of 


TAWS gave different values as seen in Table J.1. The 
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calculated temperature difference parameters were reduced 


Significantly compared to WinEOTDA outputs. 












WIinEOTDA DETECTION DELTA T (K) TAWSSeEBLECTION DELTA T (K) 
(4 km transmissivity = < a (243 Timeransmissivity = 0.589 
absolute humidity = acsclueernumiaity = 9.01) 


(90 deg) (0 deg) (90 deg) (0 deg) 


Table J.1 ~—- WinEOTDA And TAWS Detection Delta T Outputs 
Comparison For The Same Scenario Input Parameters For 
Different Sensor Altitudes 













SENSOR 
BETGHY 
CEE?) 








However, in the case of detection ranges, NFOV 
predictions were found to be the same except for the 2000£t 
sensor altitude. As can be seen in Table J.2, =TAWS 
calculated different detection ranges for varying aspect 
angles and sensor altitudes in WFOV detection. This can be 
accepted as an improvement to the insensitivity of WinEOTDA 


to changing aspect angles and sensor altitudes in WFOV 


@etection. 






















TAWS MRT DETECTION RANGE 
(4 km transmissivity =0.589) 
(Km) 


(90 deg) (0 deg) (90 deg) (0 deg) 
mor | weov | moved weov | mrova) wrov | Nrov | arov 
| 500 | 44.1 | 18.4 | 44.1 | 18.4 [ 44.1 [| 36.1 | 43.2 [ 30.8 
eer Ce mao | eo aoe eho eee 35.6 | 42.5 [ 30.7 | 
10a | ec Scene oo oeeds6 | 55.5 | 38.6 | 


Table J.2 - WinEOTDA And TAWS MRT Detection Range 
Comparison Table With The Original Scenario Parameters For 
Different Sensor Altitudes 


WinEOTDA MRT DETECTION RANGE 
(4 km pana a = 0.60) 





SENSOR 
AE IGHT 
oe) 
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